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j^JT)' Abstract 

A general canonical formalism for discrete systems is developed which can handle vary- 
ing phase space dimensions and constraints. The central ingredient is Hamilton's principle 
function which generates canonical time evolution and ensures that the canonical formalism 
reproduces the dynamics of the covariant formulation following directly from the action. We 
' apply this formalism to simplicial gravity and (Euclidean) Rcgge calculus, in particular. A 

O" 1 discrete forward/backward evolution is realized by gluing/removing single simplices step by 

5^ ' step to/from a bulk triangulation and amounts to Pachner moves in the triangulated hyper- 

surfaces. As a result, the hypersurfaces evolve in a discrete 'multi-fingered' time through the 
full Regge solution. Pachner moves are an elementary and ergodic class of homeomorphisms 
and generically change the number of variables, but can be implemented as canonical trans- 
J> ' formations on naturally extended phase spaces. Some moves introduce a priori free data 

which, however, may become fixed a posteriori by constraints arising in subsequent moves. 
The end result is a general and fully consistent formulation of canonical Regge calculus, 
thereby removing a longstanding obstacle in connecting covariant simplicial gravity mod- 
els to canonical frameworks. The present scheme is, therefore, interesting in view of many 
approaches to quantum gravity, but may also prove useful for numerical implementations. 
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1 Introduction 
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^ ■ Over the past decades, discrete gravity theories and models have been employed as useful tools 

in the study of numerous aspects of both classical and quantum gravity. The discrete structures 
usually provide, on the one hand, formulations which are amenable to numerical investigations 
[TJ[2] and, on the other hand, a regulator which is particularly convenient for the construction and 
definition of (UV-finite) quantum gravity models [3j HI [5j El [9] . An auxiliary discretization 
of a continuum theory can, however, introduce a number of discretization artifacts as well as 
break continuum symmetries, in particular, the continuum diffeomorphism symmetry of general 
relativity is generically broken in discrete gravity [3j HI [TO], [TTJ [121 EB E] • 

The vast majority of research in discrete gravity has been performed in the covariant setting, 
i.e. by means of formulations directly following from an action. Specifically, in the quantum 
theory one has thereby benefited from the many numerical methods devised for lattice gauge 
theories which are based on Euclidean path-integral techniques [HHHH]- ^ n con t ras t to this, there 
have been few efforts to construct canonical formulations of discrete gravity; first attempts |15j 
were based on a discretized space, but a continuous time. However, consistency requires that 
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the canonical theory reproduces exactly the dynamics and (possibly broken) symmetries of the 
covariant formulation, i.e. solutions of the canonical theory must replicate those following directly 
from the action. In particular, a canonical framework consistent with the discrete action ought 
to yield a discrete time evolution. A first consistent canonical formulation of simplicial gravitjQ 
has only appeared recently [101 [TT] and was based on ideas from the 'consistent discretizations 
program' [16] and an evolution scheme presented in [T7j. However, this canonical evolution 
scheme was only applicable to a special class of triangulations with both fixed topology and 
connectivity of the 'spatial' triangulated hypersurfaces and no consistent canonical formulation 
of simplicial gravity has been constructed in full generality thus far. 

One of the major obstacles in the quest of a general canonical formulation of simplicial 
gravity has been the issue of foliating a D-dimensional covariant solution by (D — l)-dimensional 
triangulated hypersurfaces which are generically comprised of different numbers of subsimplices 
which, in turn, carry the variables of the theory. The dimension of configuration or phase spaces 
associated to such discrete slices must therefore vary and a canonical evolution scheme thereby 
calls for mappings between phase spaces of different dimension. 

Why should one attempt to construct a general formulation of canonical simplicial grav- 
ity in the first place? The general answer is that classically this provides a better notion of 
discrete dynamics and time evolution, and, furthermore, in quantum gravity one is interested 
in, e.g. transition amplitudes between different three-geometries which generally requires to 
be able to identify the corresponding quantum states as elements of a Hilbert space of three- 
geometries. But there are further practical reasons why it appears useful to be able to treat a 
changing/evolving lattice with varying numbers of physical and gauge degrees of freedom in a 
canonical framework. For instance: (i) There are plenty of physical situations whose description 
may be facilitated via an adaptation of the discretization in time. For (discrete) gravity the 
prime example is an expanding or contracting universe, while a more extreme example is the 'no 
boundary' proposal |18j . (ii) Such a framework may foster the scope of numerical implementa- 
tions, (iii) It may be advantageous for quantization and help in comparing or even connecting 
various approaches to quantum gravity. Most of the latter can be roughly divided into two 
categories, namely canonical continuum quantizations (e.g. Loop Quantum Gravity |19[ I20| ) 
and covariant path integral formulations, regularized by a discrete lattice structure (e.g. Causal 
Dynamical Triangulations . Spin Foams {8j[9l 20], Quantum Regge Calculus [3]...). A link be- 
tween such approaches remains elusive thus faro It seems that a consistent canonical framework 
for simplicial gravity is a prerequisite for at least comparing these two categories, especially in 
view of Loop Quantum Gravity, on the one hand, which involves changing 'spatial' graphs, and 
Spin Foams, on the other hand, which are based on triangulations. 

Motivated by these potential benefits, it is the goal of the present article to overcome the 
technical obstacles first at the classical level by devising a canonical (discrete) evolution scheme 
for triangulations which 

1. is general and applicable to arbitrary triangulations (of fixed 'spatial' topology, but varying 
'spatial' hypersurfaces) satisfying the equations of motion following from an action, 

2. is local (i.e., for practical reasons, one may restrict oneself to a finite number of equations), 
and 

3. may be interpreted entirely from the perspective of the (D — l)-dimensional hypersurfaces. 



1 By 'simplicial gravity' we mean a gravitational theory based on simplices and triangulations, as opposed to 
the more general term 'discrete gravity' which may encompass also other discretization schemes. 

2 With the exception of the Loop-Quantum-Gravity- and Spin- Foam-formulations of the topological 3D-theory 
which were shown to be equivalent [21] . 
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The underlying basic idea is to build up a .D-dimensional triangulation (satisfying the equations 
of motion) step by step by gluing at every discrete evolution step one single D-simplex onto the 
{D — l)-dimensional hypersurface of the previous piece of triangulation. Backward evolution 
is analogously obtained by removing single D-simplices from the bulk triangulation. In such a 
way — in analogy to the evolution of hypersurfaces in canonical general relativity — a 'spatial' 
triangulated hypersurface evolves in a discrete 'multi-fingered' or 'bubble' time through the full 
D-dimensional solution. As we shall see, the gluing/removal moves which we permit can be 
interpreted within the hypersurfaces as Pachner moves, which are an elementary and ergodic 
class of topology preserving moves, mapping between any triangulations of the same topology 
by finite sequences. This local evolution scheme has to be constructed consistently in such a 
way that the equations of motion are satisfied at every step and requires the discrete action to 
be additive. The latter is fulfilled for a particular discretization of general relativity, namely 
length Regge Calculus, for which we will develop this evolution scheme here in 3D and 4D. 
However, the scheme is more general and, through a suitable change of variables, applicable to 
other simplicial gravity theories (with additive action) as well. 

In order to incorporate these ideas in a canonical framework, we must first of all extend 
the theory of discrete dynamics from the regular [22] (i.e., for systems without constraints) 
to the irregular case where constraints necessarily occur (see also |16|). before applying this 
to simplicial gravity. The central idea of our canonical framework is to resort to Hamilton's 
principle function as a generating function for a canonical time evolution; in our irregular case, 
this time evolution map will only be defined on and between constraint surfaces. Additionally, 
equations of motion transform naturally into certain canonical constraints which permit to 
suitably extend phase spaces, associated to hypersurfaces in such a way as to implement the 
Pachner moves as canonical transformations. Because of the a priori (i.e. prior to extension) 
varying phase space dimensions, the space of initial data corresponding to a given hypersurface 
cannot, in general, correspond to the space of solutions (modulo gauge) and we generically have 
to expect a high degree of non-uniqueness. Indeed, as we shall discuss in more detail in the 
main body of this paper, not every Pachner move will invoke equations of motion and rather 
allow for data that can be a priori freely chosen at that step, but may become fixed a posteriori 
by additional constraints arising in subsequent moves. Nevertheless, the end result is a general 
and fully consistent formulation of canonical Regge Calculus. 

This new canonical framework can also be viewed as a method to generate solutions and is, 
furthermore, susceptible to quantization. In particular, it admits the advantage that, by using 
the action as a generating function, the eventual canonical quantum theory should directly 
correspond to the path-integral formulation based on the same discrete action. 

The rest of the article is organized as follows: section [2] provides the basic ideas to deal 
with systems with changing phase space dimensions and summarizes all necessary facts for 
generating functions of relevant canonical transformations. Subsequently, in section[3]we outline 
the elementary evolution scheme for triangulations which amounts to Pachner moves in the 
hypersurfaces and review essential background knowledge of Regge Calculus. Section [4] discusses 
in detail the Regge dynamics and implementation of Pachner moves in 3D, elucidates how to 
start from a zero-dimensional phase space and exhibits in an example the reproduction of the 
so-called tent moves in 3D by sequences of the Pachner moves, while section [5] investigates the 
same points in 4D. In order to enhance readability of the article, we decided to move more 
technical aspects to the end of the manuscript where, after an introduction to the Lagrangian 
and canonical formalism for regular discrete systems in section 16. 1| we develop in section 16.21 a 
general formalism for singular discrete systems, which, in particular, is also applicable to systems 
with varying phase space dimension. We discuss the relevant symplectic structure and specify in 
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which sense symplectic forms are preserved, before applying this formalism to Pachner moves in 
section [631 Section [7] offers a brief interpretation of the constraints arising in the new framework 
and, finally, section [8] contains the conclusions and an outlook to further work. 

The application of the present framework to the 4D linearized theory will appear in a com- 
panion paper [23] . 

2 Canonical discrete dynamics 

The central idea |10[ \W[ \22\ [24"1 [25j [26] for generating a canonical discrete dynamics is to employ 
Hamilton's principal function S as a generating function for a canonical transformation, which 
implements time evolution. 

Hamilton's principal function S is the action S evaluated on solutions. That is, we assume 
that the action defines a well defined boundary problem, with boundary data given by some 
configuration variables Xi n i and Xfi n associated to an initial and final boundary. (We will 
assume that we can split the boundary in two parts.) Hamilton's principal function S is thus a 
function of these boundary data X{ n i and Xfi n . Since it arises by integrating out (i.e. solving for) 
all the bulk variables, we might use the term 'effective action', instead of Hamilton's principal 
function. 

Furthermore, we assume that the action is additive in an appropriate sense. More precisely, 
the action associated to a region which is comprised of two regions should be the sum of the 
actions associated to each of the two regions 

S({x} AuB ) = S({x} A ) + S({x} B ) , (2.1) 

where {x}c denotes all the dynamical variables associated to the region C . As a consequence 
of the fact that actions normally arise as space-time integrals over Lagrangians, additivity can 
usually be obtained by taking into account boundary terms. 

If the regions are such that the boundary variables of A, B are {xj n j , sinter} and \x{ n t er , x fm} i 
respectively, we find the following convolution property of Hamilton's principal function 

S{xi n i,Xfi n ) = extr Xinter 

Here 'extr' indicates that we look for the value of Xi n t er that solves the variational problem of 
the functional on square brackets in t|2.2j) . i.e. a solution to the dynamical problem defined by 
the action 

3- 'inii^- 'inter iX j 'in) • — "5(^mi) Winter) ~i~ S(Xi n t er 5 2</in) • (^-3) 

The property (|2.2p can be proven by splitting the variational problem over all variables associated 
to the region AL)B (but keeping Xi n i,xji n fixed) into three parts: one involving only the variables 
inside A, the other one only involving the variables inside B and, finally, varying with respect 
to the boundary data Xi n t er - 

The action on the left hand side of (|2.3p can, in fact, be understood as a discrete action of 
a problem with two time steps. Usually, we do not have Hamilton's principal function available 
(as this requires the solution of the continuum problem), however, one can approximate it for 
small time intervals with some discrete (one time step) action, which is a function of initial and 
final configuration data. Enumerating the time steps with a discrete label k, we will denote such 
a choice by Sk '■= S(xk-i, Xk)- 



S\Xinii Xi n i er ) + S(Xi n t er , Xfi n ) 



(2.2) 
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In order to define a Legendre transform and time evolution from such a discrete action, 
we can appeal to the standard formalism for Hamilton's principal function. Hamilton's prin- 
cipal function is also a generating function of the first kind (i.e. depends on the old and new 
configuration coordinates) and thereby determines the canonical time evolution: 

-J- 1 , V:=|^ • (2.4) 

ox k ^i dx k 

(Recall that Sk := S(xk-i,Xk)-) Here we introduced forward ~p and backward + p momenta, 
because in a discrete formulation (|2.4p allows to define two momenta at every time step: 

dx k ' dx k 

Note that the requirement + p k = ~p k , which we term 'momentum matching' implements the 
equations of motion (in the sense of (|2.2p ) for the variables Xf~. Or, conversely, the equations of 
motion implement a momentum matching of forward and backward momenta such that there 
are unique momenta for the variables at step k. Henceforth, we will often omit the superindices 
+ and — at the momenta, implicitly assuming that momentum matching holds. 

Accordingly, we can define two Legendre transformations from Qk-i X Qk to the phase spaces 
Vk-i and Vk |22j . where Qk is the configuration space at step k, 

¥ + S k : (xjfc_i,Xfc) i-> (x k , + p k ) = ( x fc>^~ 

F~S k : (x k -i,x k ) i-> (x k -i, ~p k ~ l ) = \Xk-U ~^~) ' 

These Legendre transformations are, in general, not isomorphisms because in the discrete we 
will allow for the general situation dim Qk-i ^ dim Qk- Constraints in phase space arise if the 
image of the Legendre transform is not given by the entire phase space, but by a submanifold, 
which is called the constraint manifold. Notice that here we have two Legendre transformations, 
consequently, we will encounter two different kinds of constraints, which we will call pre- and 
post-constraints. Pre-constraints result from ¥~Sk and constitute the conditions that the time 
evolution step from (k — 1) to k can take place. (In other words, that momentum matching 
can be applied at time step (k — 1).) Post-constraints, on the other hand, arise from ¥ + Sk and 
are relations on phase space that are automatically satisfied by the momenta + p k (for all initial 
values) after having performed an evolution step from (k — 1) to k. By momentum matching 
these post-constraints will provide conditions for the momenta ~p k , that is, the momenta for 
the next evolution step. The symplectic structure in the presence of these constraints will be 
discussed in detail in section El 

The formalism outlined here can be easily applied to a discrete time evolution; at least to 
the case where space-time can be foliated into disjoint hypersurfaces, which, moreover, carry the 
same number of variables. However, for a triangulation this need not be the case. Furthermore, 
we are striving for a local notion of time evolution, i.e. one in which only a small region of a 
hypersurface is evolved, or pushed forward in time. Hypersurfaces will, consequently, in general 
overlap. The goal to devise a general discrete time evolution scheme therefore gives us the 
following problems to tackle: 

(a) A time step may involve bulk variables. This will not severely complicate the situation, 
since the Hamilton-Jacobi formalism, as described here, automatically takes care of this. 
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(b) Different time steps, i.e. hyper surf aces defining instants of time, may involve coinciding 
subsets of variables. In other words, two hypersurface with different time labels, may 
partially overlap. This is an example of the 'multi-fingered' time evolution encountered in 
diffeomorphism invariant theories. This issue can also be handled easily; here the additivity 
of the action plays an essential role. 

(c) The number of variables may actually differ from time step to time step. As a result, the 
phase space dimensions change in general and we have to reconsider 'canonical evolution'. 

Let us discuss these issues point for point. 

(a): Assume we can split the configuration variables Xk at step k into two sets: 'true boundary 
variables' x\ and 'internal (bulk) variables' x\. Such a situation can be easily constructed by 
starting from some elementary evolution scheme with time steps labelled by k € % and then 
defining an effective evolution scheme where we consider only every ra-th time step (or consider 
time steps of different basic lengths) and henceforth label these rearranged steps by a different 
label n € TL. That is, for these steps we then have x l n and x l n . For example, figure Q] depicts 
a triangulation which has been built up by simplices which can be grouped up into fat slices 
bounded by non-intersecting hypersurfaces E n . Individual simplices are enumerated by A: € Z, 
while fat slices are enumerated by the new label n € TL. Now assume this to be a Regge 
triangulation where the lengths of edges are the configuration variables and assign lengths of 
edges inside fat slice n and of edges contained in X] n to step n. Thus, at step n we have the 'bulk 
lengths' l % n of edges which reside inside the fat slice n, as well as the 'true boundary lengths' 
of edges e C S n . Consequently, the action associated to fat slice n + 1 does not depend on l % n 
and we have the dependence S n+ i{l^, l l n <i, etc. 

, E n _|_i 




Figure 1: Schematic illustration of a fat slicing of a triangulation. The triangulation can be built up 
step by step by single simplices where we count such elementary steps by k E TL. The elementary steps 
can be grouped up into fat slices which we now count by n G TL. 



In such a situation, we can generally require that the action S n (x n ~i,x n ) does not depend 
on the internal variables x l n _i, as these variables only appear between times (n — 2) and (n — 1). 
The time evolution equations are then 

_ n -l ._ dff n + dSn 

Pt ' ' Pt ' dxi 

Pt dxU ~ ' P * - 9< • (2 ' 7) 

Hence, we have ~p™~ 1 = and, by momentum matching, also = (which may lead 

to further 'secondary' constraints). If the same kind of splitting into internal and boundary 
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variables occurs for all time steps re, we will obtain pf = for all n. The equations of motion 
for the internal variables are then implemented via momentum matching 

° = ~ P *= +P *=M ■ (2 - 8) 

This situation will always arise if some class of variables only appears in the action S n for only 
one label re. 

On the other hand, here we see the first occurrence of a pre-constraint; the constraint 
hypersurface is defined by pi = 0. This is a general feature: constraints appear as equations of 
motion which involve only canonical data from one time step. By redefining the time label of 
variables (here by summarizing several time steps enumerated by k into one labeled by n), one 
can transform equations of motion into constraints. 

The internal variables can be integrated out and we can define an effective action S n which 
only depends on the 'true boundary' variables x t n _ 1 and x l n . One can easily show that this leads 
to an equivalent time evolution for the remaining variables. 

(b): Assume we have the same variables x b appearing in different time steps, for instance 
x k+l = x b k for some set of variables x b (we return here to a counting of time steps by k). This 
generally happens if we implement a local time evolution, i.e. if the hypersurface in every time 
step only changes in a specific region. 

First let us consider the case that x b is not dynamically involved in the time evolution at all, 
that is, we want to implement the evolution equations 

4+1 = 4 , P k b +1 = Pb ■ ( 2 -9) 

It is not possible to implement these equation by using the action as a generating function of the 
first kind, as the fact that the x b are not dynamically involved means that neither x k nor x k , 1 
appear in the action S k+ i. However, for these variables we can use the identity transformation 
generated either by a generating function of the second (depending on old configuration and new 
momentum variables) or the third kind (depending on new configuration and old momentum 
variables): 

G 2 (xU +1 ) = ~4p k b +1 , P k b =-^i=P k+1 , x^=-^% = xt 

oxl °P\, 

rVr fc+1 n k \ - T b r> fc 7i fc+1 - - n k T k - 9Gs - T b O 1 (\\ 

<^3{ x b 'Pb) — x k+lPb i Pb — ^Tb — — Pb i x b — -jTY — X k+1-K z - W ) 

° x k+l °Pk 

Another case which will appear is that some configuration variables do not evolve x e k = x k+1 , 
while the associated momenta, however, change as either x e k or x e k+l appear in the action Sk+i 
(but not both for the same index e). That is, in accordance with (|2.6p . we want to implement 
either 

n k _ _fc+l _ 9S k+1 (x k ) k+1 _ k dS k+1 {x k+1 ) 

Pe - Pe rs e OI p e Pe ' p.* ' \ L - LL ) 

UX k UX k+l 

which we call 'momentum updating'. This can still be implemented via a generating function 
of the second or the third kind by simply adding either G 2 or G3, respectively, to the actions 
Sfe+i- As Sk+i either only depends on the old configuration variables x k or only on the new 
configuration variables x k+1 , the type of the generating function, in fact, does not change. One 
can check that these generating functions lead to the evolution equations we wished for. 
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Note that with these definitions, all momenta at any given time step can also be defined by 



V 



dS k -(x k ) 
dxh 



+ p k ds k+{xi 



dxi 



(2.12) 



where S k ~{x k ) denotes the action associated to the region which lies to the future of the hy- 
persurface with label k and S k+ (x k ) the action associated to the region in the past of this 
hypersurface (see figure [2]) . Here x k denotes all variables associated to and the future re- 
gion, while x~l denotes all_yariables associated to and the past of it. Moreover, x k can be 
any 



configuration variable]! 



Future 



Sk- 



Figure 2: Hypersurface Sfc separating 'past' and 'future' region at step k. 



This will also include the case where, say, variable x appears in the past of time k, but 
not in the future. The equation of motion is then automatically implemented by ~p k = and 
momentum matching. In other words, in general we can actually add variables which are not 
'true dynamical variables', i.e. which are not associated to the hypersurface but rather to the 
future or past of it. Then one of the associated momenta + p k or ~p k will be zero, thus enforcing 
the equations of motion by momentum matching + p k = ~p k . 

This allows us to deal with problem 

(c), namely the case where the number of configuration variables and, therefore, the dimensions 
of the phase spaces associated to different hypersurfaces Sfc-i and differ. We can formally 
extend the phase spaces by including any variables into the phase space at time (k — 1) which 
only appear at time k but do not have a natural associated variable at time (k — 1), and vice 
versa. Then we can have phase spaces of equal dimensions at the two time steps and we can 
implement time evolution by a standard canonical transformation. 

Say, we have a new variable x k+l which appears at time step (k + 1) but not at time step k 
(for a similar situation in 3D Regge Calculus see figure [3]) and the action Sfc+i only depends on 
x k+i but not on x k . Accordingly, we extend the phase space at time k by the pair (x^,p^) and 
may use the action as a generating function of the first kind, resulting in the evolution equations 

_ k _ dS k+1 _ + k+1 _ dS k+1 



3 In (|2.12[) . if the same variable x appears with multiple time labels = Xk+i, ■ ■ ., it is understood that the 
action is expressed only as a function of x = x k . 
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The variable x^ remains undetermined, for it does not appear at all in the action Sk+i- 
We may fix its value to = i^m, such that it appears as an initial datum for time step k, 
determining the data at time step k + 1. Furthermore, in (|2.13p . the pre-constraint ~p^ = 
shows up as a constraint on the variable p!^, which appears at time k for the 'first time'. Hence, 
by momentum matching it does not place any restrictions on the dynamical variables at step k. 
(If one also extended the phase spaces at previous time steps by this variable, the momentum 
p n would also be vanishing.) Additionally, we will also encounter a post-constraint: Assuming 
that N further variable pairs (x\.,p^) are involved, the time evolution will map from a (2N + 2)— 
dimensional phase space to another (2N + 2) -dimensional phase space (using that x% = x^ +l ). 
However, this map is not defined on the full phase space, because the pre-constraint p\ = 
has to hold. As a consequence, the image of the map can be maximally (2iV + 1) -dimensional, 
implying the occurrence of also a post-constraint. Indeed, in the examples in the sequel it will 
turn out that the second equation in ()2. 13[) only involves variables from time step (k + 1) and, 
hence, constitutes a constraint. 

Conversely, the case of an old variable x° k that does not have an equivalent at time [k + 1) 
can be treated analogously, that is, we extend the phase space at time (k + 1) by the pair 
( x k+i> Po +l ) ■ This time we obtain a post-constraint = 0. Likewise, we will also find a 

pre-constraint, which in the case of the Pachner moves will always be given by 

-rf - -ff • P-") 

i.e. this equation will only involve variables from time step k. Via momentum matching, ~p^ = 
+ Po, this constraint effectuates the equation of motion dS/dx° = 0, as x° appears in Sk+i for 
the 'last time'. 

In the following section, we will implement Pachner moves as time evolution maps in sim- 
plicial gravity and, with the help of these examples, explain the concepts encountered here in 
greater depth. 



3 Pachner moves as time evolution maps 

In the previous section, we presented a general method to derive a canonical time evolution 
scheme from a given discrete action. The advantage of this method is that the canonical evolution 
equations will exactly reproduce the equations of motion of the (discrete) covariant formalism. 
This is, in general, different from the procedure of first deriving a continuum canonical formalism 
from the continuum action and subsequently discretizing this continuum canonical formalism 
|15j . In the case of general relativity, where diffeomorphism symmetry plays a fundamental role, 
the latter method usually leads to (classically) anomalous constraint algebras [U \12\ fl5] . 

A canonical evolution scheme for Regge calculus which reproduces the covariant equations of 
motion has been presented in [10 Hm l2"3]. However, this scheme was only applicable to a special 
class of triangulations, namely those which can be evolved from a hypersurface by so-called 
tent moves. These tent moves [T7] are special local evolution moves that do not change the 
topology and the triangulation (that is the connectivity of the triangulation) of the 'equal time' 
hypersurfaces. Hence, in this case one does not encounter phase spaces of different dimensions. 

In the present work, we will advance and devise a general canonical framework applicable to 
arbitrary triangulations (of fixed 'spatial' topology) by implementing an elementary and ergodic 
class of local evolution moves, the so-called Pachner moves [27] . Pachner moves are elementary 
in that they involve only a fixed number of simplices during each move and are, furthermore, 
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topology-preserving and ergodic, i.e. can map between any triangulations of the same topology 
by finite sequences. In particular, the tent moves (see sections 14.21 and 15. 5p . which can involve 
an arbitrary number of simplices, can, in fact, be decomposed into Pachner moves. 

More precisely, the evolution moves to be implemented below can be interpreted entirely 
within the 'spatial' hypersurfaces (as local changes) and appear as Pachner moves therein: If 
we take the (D — 1) -dimensional hypersurface X as the boundary of a D-dimensional bulk 
triangulation, the Pachner moves in £ arise by gluing a single .D-simplex to, or removing one 
D-simplex from the bulk triangulation [27] . Notice that here we only allow for gluing processes 
of top-dimensional simplices which identify faces of one dimension less of a given simplex with 
faces of equal dimension in the hypersurface (and analogously for removal procedures). For 
instance, figure [3] depicts the 3D example of the 1-3 Pachner move in the 2D hypersurface. All 
other Pachner moves in (D — 1) dimensions can be similarly produced by gluings or removals of 
single D-simplices in the D-dimensional bulk triangulation. 





Figure 3: 3D example: gluing a single tetrahedron onto a single triangle in the 2D boundary hypersurface 
of a 3D bulk triangulation. From the perspective of the 2D hypersurface this gluing move appears as a 
subdivision of the triangle t. That is, the move appears as a 1-3 Pachner move in the hypersurface. 
This is a specific example of the situation described earlier where the numbers of variables associated to 
the two hypersurfaces before and after the new move differ due to the three new edges. 



This elementary procedure provides a more compelling connection between the covariant and 
canonical picture: the discrete evolution of the hypersurface can be reinterpreted as building 
up the bulk triangulation, and, hence, the (discrete) space-time, step by step by simplices. In 
analogy to the situation in canonical general relativity, the triangulated hypersurface evolves 
in a discrete 'multi-fingered' (or 'bubble') time through the full (discrete) space-time solution. 
In this light, every gluing move at a given evolution step k which adds a simplex to the 'past 
triangulation' can also be viewed as a removal move which subtracts a simplex from the 'fu- 
ture triangulation', and vice versa (see figure [3] for a schematic illustration). Formally, this is 
guaranteed by momentum matching. 

Since Pachner moves have inverses and, thus, in order to establish a genuine implementation 
of Pachner moves into the mathematical formalism, we have to ensure that every phase space map 
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'Future' 




Figure 4: Removal moves for the 'past triangulation' are equivalent to gluing moves for the 'future 
triangulation' and vice versa. 



Pachner j corr esponding i Q the action of a Pachner move on a hypersurface has an inverse which 
corresponds to the action of its inverse move (of course, via the gluing/removal correspondence 
applied to the same simplex), i.e. Pachner f o Pachner j _ This, indeed, will be the case. 

As Pachner moves do not change the topology of a triangulation, a (D — l) hypersurface with 
topology T evolved by Pachner moves will lead to a D-dimensional triangulation of topology 
[0, 1] x T. Hence, just as usual canonical evolution schemes, also the present scheme assumes 
or implements a non-changing topology of the 'equal time' triangulations. This might be ad- 
vantageous, for instance, for quantization, if one wants to suppress topology changes of spatial 
hypersurfaces (interpreted as the production of baby universes). Specifically, the approach of 
Causal Dynamical Triangulations (CDT) [6] shows that this can lead to a much more regu- 
lar large scale limit than in the case of (Euclidean) Dynamical Triangulations [28], where 
arbitrary topology changes are allowed forQ 

On the other hand, a Pachner move evolution scheme allows for a change of the number 
of variables 'in time'. There are numerous physical situations, where it might be interesting to 
adapt the discretization (density) in time, for instance, for an expanding or contracting universe. 
In fact, we will show that, akin to the 'no-boundary' proposal [18], we can start with an 'empty 
triangulation', obtain the boundary of a L)-simplex in the next time step and evolve this to a 
bigger and bigger triangulation. 

Furthermore, Pachner moves implement the idea of a 'fluctuating lattice' in a controlled 
manner, that is, the idea that the lattice, or discretization, is not fixed in time (or space-time) 
but is either summed over or is determined by dynamical considerations. Indeed, we will see that 
even the classical dynamics may prefer or suppress certain Pachner moves and, hence, determine 
the evolution of the connectivity of the triangulation itself. 



3.1 Regge calculus 

In this article, we will develop this evolution scheme for a particular discretization of general 
relativity, namely Regge calculus [21 [29]. However, the basic ideas are general and straightfor- 

4 The evolution moves implemented here yield triangulations akin to the triangulations employed in CDT. The 
splitting off of baby universes, aka spatial topology changes, on the other hand, would require the implementation 
of moves distinct from the Pachner moves. For instance, in D dimensions one could allow for a move which identifies 
all D + 1 [D — l)-faces of a D-simplex with D + 1 (D — l)-simplices in a (D — l)-dimensional hypersurface. 
Furthermore, by only allowing top-dimensional simplices to be glued onto faces of one dimension less, we are 
disallowing singular moves which, e.g., only identify a single vertex of the new simplex with a single vertex of the 
hypersurface. 
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wardly adaptable to other discretization schemes fulfilling the basic prerequisite of additivity of 
the action; whenever we glue or remove an elementary building block, i.e. a simplex, we have to 
add or subtract the corresponding piece of action. 

To begin with, we will review some basic notions of Regge calculus based on length vari- 
ables, before studying the details of the individual Pachner moves in sections [Hand [5] below. For 
formulations using other geometric variables see, for instance, [30^ 131], Standard Regge calculus 
is based on a fixed triangulation where the basic variables are simply the length variables asso- 
ciated to the edges of this triangulation. These length variables (assuming generalized triangle 
inequalities are satisfied, which we will always do) specify a piecewise linear geometry for the 
triangulation in question. The equations of motion are determined by a variational principle 
from an action, the Regge action, which can be viewed as a discretization of the Einstein-Hilbert 
action. 

The (Euclidean) Regge action (without a cosmological constant term) for a L>-dimensional 
triangulation T with boundary dT and interior T° := T\dT is given [2j [32] by 

S = ^V h e h + V ^h , (3.1) 

hcT° hCdT 

where h stand for the (D — 2)-dimensional subsimplices (sometimes called "hinges" ) of T, that 
is edges e in 3D and triangles t in 4D. denotes the volume of the hinge h. The deficit angles 
€h and exterior angles ifth are given by 

e h = 2tt - °h for h C T° , (3.2) 

aCh 

= k h n - Y,°h ioihcOT, 

aCh 

and in both cases the sum ranges over all .D-dimensional simplices a which contain h, and 0% is 
the interior dihedral angle in the simplex a between the two (D — l)-dimensional subsimplices 
that meet at h C o". The number is fixed according to how many pieces I of triangulations 
are glued along the hinge in question: Because of the boundary term in (I3.ip . the action is 
additive if we glue two pieces of triangulations together and choose = 1. In the canonical 
evolution scheme we will glue top-dimensional simplices to bulk triangulations. That is, the 
smallest triangulations we consider is just one tetrahedron r in 3D or one 4-simplex a in 4D, 
of which the actions are just boundary terms, 

S T = ^l e (k e ir-ei) 

eCr 

S a = J2A t (ktK-6?) , (3.3) 

tCcr 

where l e denotes the length of the edge e and At the area of the triangle t. Again, k e and kt are 
to be determined by the gluing process. In the canonical evolution scheme presented below, we 
will either have k e , kt = 1 or k e , kt = 0. 

The equations of motion are obtained by varying the action (|3.ip with respect to the lengths 
l e of the edges in T° (while keeping the edge lengths in the boundary dT fixed). The boundary 
term in ([3.1 p ensures that, for all edges e G T°, the Regge equations read 

£f»-°- (^) 
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Note that in 3D the equations of motions enforce e e = 0, that is, vanishing deficit angles and 
therefore vanishing curvature. 

Finally, the variation of the deficit angles appearing in the action vanishes because of the 
Schlafli identities 

Ti^ - o 

eCr e 

J>f _ . (3.5) 

tC<7 



4 Pachner moves for 3D Regge calculus 

In the next subsections we will spell out the details of the dynamics of Pachner moves for 3D 
Regge gravity (without a cosmological constant). Before delving into the details, we will make a 
few remarks on the dynamics of 3D gravity. As mentioned in the previous section, the equations 
of motion for gravity require that all deficit angles vanish and thereby that the piecewise flat 
geometry is flat. In the canonical description we describe the dynamics of the 2D hypersurface, 
that is, the changes in the intrinsic and extrinsic geometry of the hypersurface evolving through 
the 3D space-time. As the latter is flat, we will essentially describe different embeddings of a 
2D hypersurface into flat 3D space. 

This allows us to determine Hamilton's principal function: since the deficit angles vanish for 
solutions, only the boundary term remains 

s(n = J2 ieMn ■ (4.i) 

eCdT 

Here the l e are the lengths of the edges in the boundary of the triangulation. The momenta 
are given by the derivatives of Hamilton's principle function with respect to the edge lengths. 
Because of the Schlafli identity, the terms with derivatives of the exterior angles vanish and 
we obtain (the different signs in (|2.4j) can be taken into account by taking the appropriate 
orientation of the exterior angles) 

Pe = ^e ■ (4.2) 

Thus, the momenta conjugate to the length variables are given by the exterior angles. The 
same result was obtained in [33] starting from a first order (connection) formulation of discrete 
gravity. 

As a consequence of the fact that the Pachner move dynamics allows to add degrees of 
freedom during the evolution, we can start from a small triangulated hypersurface and evolve 
to a much bigger one. In particular, we can start from an empty triangulation at time k = 0. 
This can be seen as just a technical method, however, might also be useful (in 4D) to explore 
the 'no-boundary proposal' [15] . 

In order to produce a triangulation with a hypersurface of, e.g., spherical topology, we can 
start with an evolution move from the empty triangulation to the boundary of a tetrahedron! 
A tetrahedron has six boundary edges Z", n = 1, . . . , 6, hence, the phase space at time k = 1 
is 12-dimensional. Accordingly, we can extend the phase space at time k = by the six pairs 

(«)• 

5 We could equally well begin by producing initial hypersurfaces of more complicated topology. 
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Using the one-tetrahedron action as generating function of the first kind 

Go-i(M) = = TO)) , (4-3) 

nCr 

the equations of motion amount to 

P°n = , pi = = vr - B (JJ) • (4.4) 

That is, the Z™ remain undetermined. However, we can set Zq = Z™ with the understanding 
that the Z" are initial data which appear only at time step k = 1. The momenta at Zc = 
are constrained to vanish. But also at time k = 1 we obtain a post-constraint because all the 
momenta are given as functions of the length variables at time k = 1. Indeed, this equation 
just expresses the fact that the momenta are exterior angles, which for a flat tetrahedron are 
determined by the intrinsic geometry of the boundary surface, i.e. the length variables. Hence, 
we have a totally constrained system, as all the momenta are determined as functions of the 
length variables. 

This hypersurface can be evolved into a more complicated hypersurface (with spherical topol- 
ogy) by means of the Pachner moves described below. Momentum updating will then take care 
of the geometrical updating of the exterior angles, i.e. the exterior angles change according to 
the dihedral angles that are added/subtracted by gluing/removing tetrahedra onto/from the 
hypersurface. 

The dynamics of 3D gravity is special because the fact that the system is totally constrained 
does not even change for more complicated triangulations (of spherical topology); through con- 
straints all momenta at time k will be determined by the length variables of step k. These 
constraints express the fact that we are always dealing with a 2D hypersurface embedded into 
3D flat space and that, furthermore, the momenta are the exterior angles. Accordingly, if we 
considered a parallel transport of a 3D vector along a small loop around a vertex v of the hyper- 
surface, we should obtain an identity transformation. This parallel transport can be expressed 
[Ml [35] as a sequence of rotations 

P v = R(a eie2 )R(ij e2 )R(a e2ea )R(^ ea ) ■ ■ ■ R(^ ei ) = Id , (4.5) 

where e±, &i, . . . denotes some cyclic ordering of the edges around the vertex v, R{ot eiei+1 ) denotes 
the rotation in the plane spanned by the two edges e^ej+i with the angle a eiBi+1 and R{ip ei ) 
denotes the rotation around the edge by an angle tfj ei . Note that the exterior angles equal the 
momenta and that the interior angles a ee / can be expressed as functions of the length variables 
(in the 2D star of the vertex). Hence, the condition that the 3D parallel transport matrix should 
be the identity gives us three (as the matrix is in 50(3)) constraints on the phase space data for 
every vertex in the hypersurface. For a triangulation of spherical topology we have 3(tu = Jje + 6 
for the number of vertices (ft; and the number of edges jje. Therefore, we have at least as many 
constraints as configuration (or momentum variables). In fact, there are six more constraints 
than edges, because there exist six relations between the constraints (as can be checked explicitly 
for the example of the tetrahedron) . These six relations correpsond to the three global rotations 
and three global translations which change the embedding of the 2D triangulation in 3D flat 
space, but do not change any of the geometrical data, i.e. neither lengths or exterior angles. 

The constraints (|4.5f) will be preserved by the Pachner moves, as these Pachner moves will 
implement the equations of motions, i.e. flatness of the 3D triangulation. Furthermore, momen- 
tum updating will ensure that the momenta are always given by the exterior angles of the 2D 
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hypersurface. Consequently, the canonical data at every time step will describe a 2D triangu- 
lation embedded into flat 3D space for which the relations (|4.5p hold. Moreover, as we shall 
see shortly, the 1-3 Pachner move generates one vertex, and in conjunction with this vertex 
also three (post-) constraints of the form (|4.4p . which are just a rewriting of the form (|4.5p 
for three-valent verticesjf] The 2-2 Pachner move and the 3-1 Pachner move will not generate 
vertices and therefore also no additional constraints. The dynamics prescribed by these moves 
will, however, preserve the constraints of the form (|4.5p for the reasons just given. 

As a mnemonic, for the description of the Pachner move dynamics in both 3D and 4D we 
will use the following edge indices in order to label and appropriately distinguish the various 
length and momentum variables: 

• e labels edges contained in the D-simplex of the Pachner move which occur in both 
and Efe+i, 

• n labels new edges introduced by a Pachner move which occur in Efc+i but not 

• o labels old edges removed by a Pachner move which occur in but not Sfc+i, 

• b labels edges contained in both and S^+i which are not involved in the Pachner move 
and from that perspective may be considered as boundary edges. 

4.1 The 1—3 Pachner move 

Consider a 3D triangulation with a boundary that we will consider as a 2D hypersurface at 
time k. Glue to this boundary a tetrahedron r such that one of the triangles is identified with 
a triangle t in the 2D hypersurface We obtain a new boundary T,k+i- From the perspective 
of the hypersurface this gluing can be interpreted as a 1-3 Pachner move, i.e. the triangle t 
is replaced by three triangles, that share one vertex v in the middle and have the same ID 
boundary, consisting of three edges as the original triangle t (see figure [3|). 

Note that the tetrahedron can be glued with two different orientations to the hypersurface. 
This can be interpreted as gluing the tetrahedron either on the upper side (with a future pointing 
tip) or on the bottom side (with a past pointing tip) to the hypersurface, or, alternatively, as 
gluing or removing a tetrahedron to or from the bulk triangulation, respectively (see figure 
[5]). The Regge actions for the different orientations of the tetrahedron just differ by a global 
sign (which can be understood to arise from the oriented exterior angles). This agrees with 
the interpretation of gluing or removing a tetrahedron with, say, positive orientation: if we 
remove a tetrahedron from the bulk triangulation, we would have to subtract the action for the 
positively oriented tetrahedron, the alternative is to add the action (glue the tetrahedron) with 
the negative orientation. 

Every Pachner move in 2D and 3D can be interpreted as gluing a simplex in one of the two 
different orientations to the 3D and 4D bulk, respectively. These two possibilities can, in general, 
also be seen as gluing and removing a simplex with, say, positive orientation, respectively. We 
will henceforth assume that the two possibilities are encoded in the orientations and therefore in 

6 The difference between the two forms of constraints is that (|4.4[) is linear in the momenta, whereas (|4.5|) 
involves cosp e - Indeed, the constraints should rather involve the square (or the cos) of the momenta as in the 
continuum. This, however, can also be taken into account in (|4.4[) : Constraints quadratic in the momenta indicate 
that the constraint hypersurface has two pieces corresponding to the two roots of the quadratic equations. The 
two pieces correspond to the possibility of allowing for both orientations of the tetrahedron (i.e. allowing both 
signs for the exterior angles). We can replace (|4.4[) by cosp„ = cosf/> n to take care of this fact. 
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a) b) 

Figure 5: a) The 1-3 gluing Pachner move, b) The 1-3 removal Pachner move. 
The clashed edges are the three new edges. 

the signs of the exterior angles in the action fl3.3|) . Thereby we can summarize these two cases 
into just one, which we will mostly refer to 'as gluing a simplex' to the bulk triangulation or 
hypersurface. 

In the 1-3 Pachner move one triangle t is replaced by the same triangle subdivided into 
three new triangles (see figure E]). Therefore, for this Pachner move we will have edges of three 




Figure 6: The 1-3 Pachner move and its inverse, the 3-1 Pachner move. 

different types: edges e& not participating in the dynamics, three edges e in the boundary of the 
triangle t, for which ^ = and the three new edges e n with lengths which only appear 
at time step (k + 1), but not at k. We will regard the action 

S T = ^^ +1 Hl(/E +1 ,^ +1 )) + ^^ +1 (vr-^(/I +1 ,^ +1 )) (4.6) 

eCT nCr 

associated to the tetrahedron glued to the hypersurface as a function of the three new edge 
lengths and the three edge lengths Note that we fixed the factors k e appearing in 

(|3.3I) to k n = 1 for the new edges l n and to k e = for the 'boundary edges' l e . This will also 
work for the other moves, we will choose k e = for all boundary edges and k n = 1 or k Q = 1 for 
edges which either appear or 'disappear' during the evolution move. In this way the ir factors 
add up correctly to 2ir after a sufficient number of evolution moves. 

According to the discussion below (12. lip , we add the generating function G^{l e k+1 ,p^) for 
the edges of type e and b to the action. Hence, the generating function is 

^1-3(4+1)^6! lt+l)Pe> ^ki^k+l) = 4+1^6 + / l ^k+lPe + ^T^jfc+i) ffi+l) • (4-7) 

b e 
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The equations of motion are then 




fc+i 



Pb 



(4.8) 



Pe + 



8S T 



Pe ~ ^e('fc+l, ifc+l) 



(4.9) 



fc+i 



Pe 




fc+1 _ 



8S T 



(4.10) 



P 




= 7T 



The momenta p e are just updated to agree with the exterior angles of the evolved spatial 
hypersurface. The momenta p\ at time step k have to vanish, since the Zj? are not dynamical 
variables at this time step. As was mentioned in the discussion at the end of section [21 we 
have to expect a post-constraint for every new edge variable l n . Indeed, equations (|4.10p are 
constraints, requiring that the new momenta p^ 1 are again given by the exterior angles of the 
new hypersurface (coinciding with the three exterior angles of the added tetrahedron) which, 
however, can be expressed as functions of the length variables Z|, ijZJL. 1 only. 

Notice that for this evolution step not only the ZJ? remain undetermined, but also the Z£ +1 
can be chosen arbitrarily (the generalized triangle inequalities have to be satisfied though). We 
will set Z£ = Z^ +1 , so that the Z£ can be interpreted as initial data that determine the data at 
time {k + 1). Later we shall see that also for some types of Pachner moves (in 4D) the edge 
lengths of new edges can be chosen arbitrarily. This freedom can be understood as a choice of 
initial data, which becomes only relevant at the time step, at which the new edges appear. On 
the other hand, we will also see (in 4D) that pre— constraints, appearing in consecutive evolution 
moves, might fix these edge lengths a posteriori. 

4.1.1 The 2—2 Pachner move 

Consider the situation in which a tetrahedron r is glued to in such a way that two of its 
triangles, t\ and t%, and thus five of its edges are identified with two neighbouring triangles 
and their five edges in This gluing move (with positive orientation) is only possible if the 
extrinsic dihedral angle ip^ at the edge e = o along which the two triangles in are identified 
is negative. If it is positive, it is only possible to remove the corresponding tetrahedron r from 
the triangulation, or, alternatively, to add a tetrahedron with negative orientation, see figure 
[7J From the perspective of the hypersurface, these elementary moves appear as 2-2 Pachner 




Figure 7: The gluing/removal 2-2 Pachner move involves two triangles t\ and fa in hypersurface £fc. 
moves. That is, two triangles sharing an edge e = o are replaced by two triangles sharing a new 
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2-2 



n 



2-2 



Figure 8: The 2-2 Pachner move is its own inverse. 



edge e = n (see figure [8]). Note that while this move removes one edge o from the hypersurface 
and, instead, introduces the new edge n, it does not introduce a new vertex. Again, the two 
pairs of triangles have the same boundary of four edges e and four vertices. The edges o, n are 
transversal, that is o and n are connecting the two opposite pairs of vertices. Additionally, the 
2-2 Pachner move is its own inverse, which from the 3D perspective is taken into account via 
the gluing/removing convention (and the different global signs in the action for the cases with 
different orientations) for the Pachner moves. 

There are four kinds of edge lengths: apart from ijg/k+v^k/k+i anc ^ one l en §th we a ^ so 
have one length l k which appears only at time step k but not at (k + 1). Thus, we will extend 
the phase spaces at time k and (k + 1) by the pairs {l k ,Pn) an d (lk+n Po +1 ) i respectively. 

Accordingly, we will choose a generating function of the first kind in the variables l n , 1° and 
of the third kind in the variables l b ,l e : 

^2-2(4+1^6; ^fc+liPe! ^J^fe+li ^fcVfc+l) = A+lPb + ^ ^k+lPe + £fc Jjfc+l) , (4-11) 

b e 

where 

+ ^fc+i ~ OnQt+lilki ^fc+i)) + 

l° k {n-9M + iX k ,l n k+1 )) ■ (4.12) 

(4.13) 

Pe-0e(l e k+ l,l%+l) , (4.14) 

Tr-e n (i e k+1 ,i° k j n k+1 ) , (4.i5) 

-vr + M^M+i)- (4-16) 

As discussed generally, we have p^ +1 = and p 1 ^ = 0. In contrast to the 1-3 move, where 
the new edge lengths l k , l k+1 were undetermined at both time step, here the edge length l k+1 
is determined by equation (|4.16p as a function of the initial data, which also involve l k at time 
k. (Again, we can define l k = l k+ i but this time the interpretation of l k as additional new data 
does not apply, rather l k is in this case constrained, that is determined by the other initial data 
at time k. The same applies to Z£ +1 if we go backwards in time.) Note that — via momentum 
matching — (|4.16p implements the Einstein equations e Q = for the edge o which becomes a 
bulk edge in the course of the 2-2 (gluing) move. Equation (|4.16p demands that ~p k Q is given by 



S r (ll +1 ,l k ,l k+l ) — (—9l(l k+ i,l k ,l k+ i)) 

eCr 



The equations of motion are given by 



P, 



l k 


— 'fc + 1 ' 


P k b +1 = 


Pb , 


l k 


— 'fc + 1 ' 


p* +l = 


k 9S T 

P " + dl* = 
0i k+l 


Pn 


= o , 


P k n +l = 


f)in 
° l fc+l 


fc+1 
o 


= o , 


Pi = 


dS T 
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the (oriented) exterior angle of the tetrahedron that is glued to the hypersurface. Momentum 
matching, on the other hand, imposes ~p\ = + Po, where + p\ is determined by the previous 
moves and given by the (oriented) exterior angle of the hypersurface determined by the bulk, 
and thus that these two exterior angles add up to a vanishing deficit angle around the edge o. 
This condition will, in general, fix the edge length as the exterior (or dihedral) angle at the 
opposite edge o of the tetrahedron depends on this length. 

The difference to the 1-3 move is that here we have both a 'new edge' and an 'old edge'. 
Hence, the argument in section [21 according to which we have to expect a pre-constraint (because 
of the 'old edge') and a post-constraint (because of the 'new edge'), does not apply, as this 
argument relied on determining the maximal dimension of the image of the time evolution map. 
Here the counting changes since we have both a new and an old variable and due to the conditions 

= and = 0, we do not necessarily expect a further constraint based on this argument 
(where we ignore that remains undetermined). A priori, with the exception of = 0, 
all momenta at time (k + 1) involve edge lengths from time k and time (k + 1), in particular, 
(335|) for P % +1 E 

The fact that there are no new post-constraints arising in the 2-2 move should, however, not 
be confused with the feature that the system under consideration is totally constrained (that 
is, for spherical topology of £ all momenta are constrained as functions of the edge lengths). 
Rather, the 2-2 move preserves the (flatness) constraints, which encode that the hypersurface in 
question bounds a flat bulk triangulation, by implementing the flatness condition for the edge 
that becomes a bulk edge. 

4.1.2 The 3 1 Pachner move 

The 3-1 move is the inverse of the 1-3 move: Consider, therefore, a three-valent vertex v in 
Efc whose adjacent three edg equipped with extrinsic dihedral angles ip which are all 

negative, i.e. the vertex is pointing into the hypersurface and represents the tip of a tetrahedron 
which is upside down. In this situation we can glue a tetrahedron r to this surface by identifying 
the three triangles sharing vertex v in r with the corresponding ones in T,^. Consider now the 
opposite situation, where all ip at the three-valent vertex are positive. In this case the vertex 
represents the tip of a tetrahedron r which is sticking out of the hypersurface and we may remove 
this tetrahedron (or equivalently glue a tetrahedron with opposite orientation) . These situations 
are depicted in figure [5] if the orientation of the arrows is reversed. From the perspective of the 
hypersurface, these elementary moves appear as a 3-1 Pachner move, see figure Prior to this 
move, all six edges involved in the move are edges in S^. During the move the vertex v, as well 
as the e Q are removed (as these become internal to the bulk triangulation). 

By virtue of the fact that there are no new edges introduced during the move, the equations 
of motion play rather the role of constraints in this case which have to be satisfied in order for 
this move to be allowed. Since up to we have always solved the equations of motion during 
the elementary evolution steps, the entire triangulation up to will be embedded in a flat 3D 
manifold. The three equations of motion of the edges which become internal, therefore, do not 
add any new constraints and are thus trivially satisfied. 

7 Since 1° remains undetermined at time (k + 1) and l n at time k, one might wonder whether it is possible to 
summarize the two variables l°,l n into just one variable (with two different time labels). Indeed, it is possible 
[36] to define a canonical time evolution in this way; the method presented here, however, is more similar to the 
1-3 and 3-1 move and, moreover, to the Pachner moves in 4D. Nevertheless, this remark shows that canonical 
time evolution maps can be redefined by just relabeling or identifying variables with each other. We will see a 
similar situation in the tent move evolution in sections 14.21 and 15.51 below, where several Pachner move steps are 
grouped together and labeled as a single tent move time step. 
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Figure 9: The 1-3 and 3-1 Pachner moves. 



The 3-1 move is the inverse to the 1-3 move, hence, we have three kinds of edge lengths: 
apart from ^/k+v^k/k+i there are three edge lengths of type 1% which appear only at time step 
k but not at (k + 1). We therefore extend the phase space associated to by three pairs 

We choose the action to be a function of the edge lengths at time k and, accordingly, use for 
edges of type e and b a generating function of second type. According to the general description, 
we have 

G 3 -i(4^V^e +1 ;M + i) = -£^ +1 -E^ +1 + S T (ltX k ) (4.17) 



with 



sMJi) = £igH£(M)) + -*o(M)) • (4-18) 



eCr oCr 

The equations of motions are then given by 



L k+1 


— l b 

— l k > 

— L k j 


p\ = 

p1 = 




= o , 


P k = 



Pt +l , (4-19) 
P" +1 ~W k =Pe +1 + e M,® , (4-20) 



dS T 



-* + e (i%,i%) . (4.2i) 



The second equation in (|4.20p defines the updated momenta p^ +l as a function of the phase 
space variables at time k, while the second equation of (|4,2ip is a pre-constraint, this time on the 
phase space variables at time k. That is, to perform this 3-1 move, this condition on the phase 
space variables has to be satisfied. Again, as for the 2-2 move these pre-constraints implement 
the equations of motions, namely flatness for the edges which become bulk edges during the 
move. In general, there are two possibilities: either i) the pre-constraints are automatically 
satisfied, if we consider a hypersurface which has been evolved in the manner described here 
by Pachner moves, or ii) the constraints are only satisfied for specific initial data, including the 
kind of additional initial data which arise, for instance, by the 1-3 Pachner move. For the 3D 
case, only the first possibility takes place as argued above. In 4D, on the other hand, we will 
generally encounter the second possibility, namely data which are a priori free to choose in a 
certain Pachner move may become fixed by pre-constraints arising in later moves. 



4.2 Example: 3D tent moves 

The Pachner moves generically change the connectivity and the number of edges in the trian- 
gulation. There are, however, combinations of Pachner moves, e.g. the so-called tent moves 
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[17\ [TT| HUj . which do not change the connectivity of the triangulations and therefore induce 
a canonical dynamics in the standard interpretation, that is, between phase spaces of equal 
dimensions. 

A tent move can be constructed by picking some vertex uJ! in a (D — l)-dimensional trian- 
gulated hypersurface S n and subsequently defining a new vertex v n +i to the 'future' of v n which 
must be connected by an edge to v n (the 'tent pole'). Denote all other vertices in S n to which 
v n is connected by 1, . . . , N. Connect v n+ i to each of these vertices 1, . . . , N by N edges. The 
new vertex i> n +i then lies in the new hypersurface and is iV-valent in S n+ i as v n is in S n 

(the tent pole is now internal) and for every (D — l)-simplex cr(v n ij . . .), i,j € 1, . . . ,N, there 
is now a (D — l)-simplex cr(v n +iij . . .). That is, the triangulations of the two hypersurfaces are 
the same. Notice that each tent move only involves the (D — 1) star of the vertex v n in E n . The 
situation for a four-valent 3D tent move is illustrated in figure [TOl 



V n +1 
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Figure 10: A 3D tent move applied to a four-valent vertex v n in a 2d Cauchy hypersurface. This four- 
valent tent move can be reproduced through a gluing sequence of four tctrahcdra which correspond to a 
sequence of one 1-3, two 2-2 and one 3-1 Pachner moves in the hypersurface. 

The new piece of D-dimensional triangulation between S n and with the N + 1 new 

edges e(v n+ ii) and e(v n v n +i) can be decomposed into a sequence of gluings of single D-simplices 
and the tent move evolution may therefore be described in terms of a sequence of Pachner moves 
in the hypersurface. In particular, in 3D simply pick an N-valent vertex in some S n and perform 
one 1-3 Pachner move, (N — 2) 2-2 Pachner moves and a final 3-1 Pachner move in order to 
generate £ n+ i. 

To define the Pachner move dynamics, we used the Regge action as (part of) a generating 
function, so that the resulting canonical dynamics leads to the same equations of motions as 
the ones obtained by varying the action. We can again employ the Regge action to define the 
canonical evolution corresponding directly to the tent moves [10\ [TT\ [23] . By construction, this 
dynamics will coincide with the one obtained by performing the sequence of Pachner moves 
(matching appropriately the edge labels). This highlights the remarks put forward earlier that 
canonical time evolution maps may be redefined by relabeling variables. For the tent move 
dynamics we can either keep the length of the tent pole as an internal variable in the sense of 
problem (a) in section [21 see equations (|2.7p . or, alternatively, integrate this length out and work 
with Hamilton's principal function as a generating function: 

To this end, notice that the 3D triangulation — the 'tent' — we are gluing to the hypersurface 

8 Note that henceforth, in order to avoid confusion, we will enumerate tent moves by n £ Z5, while the elementary 
Pachner moves into which the tent moves can be decomposed are counted by k £ 7L. 
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has two 2D boundaries given by the 2D stars of the vertices v n and v n+ \. These two 2D 
boundaries meet in a ID boundary, given by the (cyclically ordered) vertices 1,...,N and the 
edges connecting these vertices. The tent consists of N tetrahedra sharing the tent pole as an 
edge, and therefore each having vertices v n , v n+ \ and i, i + 1. The tent pole is an internal edge, 
whereas the 2D boundary has edges e(v n , i), e(v n+ i, i) and edges e(i,i + 1) with edge lengths 
l l n , l l n+ i, i = 1, • • • N and l l n +1 = l^+i , respectively. We can solve the equation of motion for the 
tent pole as a function of the boundary edges. This equation of motion will just require that 
the deficit angle around the tent pole is vanishing. Using this result in the action, we obtain 
Hamilton's principal function for the 'tent' 

C fji 7i ji,i+l Q T(v n ,v n+1 ,i,i+l) , li / Q T(v n ,v n+ i,i,i+l) Q T(v n ,v n+1 ,i-l,i)\ 

i i 

En ( ' a r(v n ,v n+ i,i,i+l) Q r(v n ,v„ + i,i-l,i)\ , . 00 n 

i 

Taking into account the additional terms Pii+i i n the generating function we obtain the 
following equations of motion 
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where again the Schlafli identity was used. Now, a priori the dihedral angles on the right 
hand side of the equations (|4.24I4.25[) depend on both sets of edge lengths l % n and l n +l . But, 
for instance, for the tent move at a three valent vertex the right hand side of (I4.25|) does not 
depend on the lengths l l n and we obtain three post-constraints. (The geometry around the 
new vertex is the one around a vertex of a flat tetrahedron and the dihedral angles can be 
expressed as functions of the six edge lengths l n +x,ln+i ■) Likewise, (|4.24p will give us three 
pre-constraints as the right hand side can be expressed as a function of l l n , only. For tent 
moves at higher valent vertices we can still conclude that we have three pre-constraints and 
three post-constraints of the form (|4.5j) . To this end one just has to use the same arguments as 
for (|4.5p applied to the triangulation of the tent, which itself is a flat triangulation. 

Hence, we will have three pre- and three post-constraints for the tent move. From this we 
can conclude that three of the edge lengths l l n+ i remain undetermined by the equations of motion 
(14.231 14.25P . (Correspondingly, there is a three parameter set of initial data at step n which can be 
evolved to the same data at time (n + 1), and vice versa.) These are the same three parameters 
that are left undetermined by the 1-3 move in section 14.11 and correspond to lapse and shift 
degrees of freedom. The three post-constraints of the 1-3 move are of the form as encountered 
here and the present discussion shows that these constraints remain preserved through all the 
additional Pachner moves which make up the tent move. Also the pre-constraints which need 
be fulfilled so that the final 3-1 move can be performed, will automatically be satisfied, if the 
triangulation has been correctly evolved by Pachner moves. 

Finally, we remark that the lapse and shift degrees of freedom which remain free in the 
tent move can be chosen to be infinitesimally small. This allows to recover a continuous time 
evolution (plus Hamiltonian and constraints generating this evolution, which coincide with (14. 5p ) 
from the discrete time evolution presented here. 
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5 Pachner moves for 4D Regge calculus 



While the 3D Regge dynamics leads to flat geometries, the 4D Regge equations allow for curved 
solutions. This renders the 4D dynamics significantly more complicated than in the 3D case, in 
particular, the preservation of the constraints will generally not hold. Related to this is the fact 
that lapse and shift degrees of freedom will, in general, not remain free but will become fixed 
by pre-constraints (unless one considers initial data which lead to flat solutions). 

The Regge action for a 4D triangulation (without cosmological constant) is given by (|3.1l3.2p 
for h = t, 

tcT° tedT 



The equations of motion for an inner edge e read 

d/. 



tDe 



which — depending on the boundary data — allow for flat solutions et = as well as solutions 
with curvature et ^ 0. 

Despite the latter, we can show that Hamilton's principal function reduces to a boundary 
term (as is the case in the continuum). To this end, multiply each of the equations (|5.2|) with 
the length l e and sum over all (inner) edges. We will employ the Euler identity (see, e.g., [37] ) 

£4f = 2At (5 - 3) 

eCt 

for the area of a triangle. Using the equations of motion (|5.2p . we can write 

- 

eCT° iCe 

= ^ £ A t e t + £ E ^ (5 - 4) 

tCT° t C T t a eCtneCT 

where T t ° denotes the set of triangles which have all three edges contained in the bulk T° and T® 
denotes the set of triangles in T° where at least one of the edges is in the boundary dT. (Note 
that an inner triangle can have either one, two or all three edges in the boundary.) Employing 
equations (I5.3j) and (I5.4j) in the action (|5.1j) . we discover that Hamilton's principal function 



eCdT 



\- dA t ^ dA t 

,tDentcT° tDelltCdT 



(5.5) 



is given by a boundary term. The form of this last expression suggests that the momenta are 
given by the combination of exterior angles and deficit angles 'near' the boundary in square 
brackets on the right hand side of (15. 5\\ . 



dA t OA 

w et + E 

Den(CT° tDentC&T 

Indeed, momentum updating during the Pachner move evolution will confirm this 
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As in 3D, we can produce a Regge triangulation with boundary topology given by the 3- 
sphere^l by considering first an evolution move from the empty hypersurface to the boundary of 
a 4-simplex. The generating function would simply be given by the one-simplex action 

G -!(M) = s a (m = , (5-7) 

t 

and gives rise to the following equations of motion 




Again, just as in the 3D case, we also obtain post-constraints here which determine all 10 
momenta at time k = 1 as a function of the edge lengths at time k = 1. Hence, the boundary 
of a simplex is a totally constrained system. This property does not change under 1-4 Pachner 
moves (yielding so-called stacked spheres) which will introduce four post-constraints for every 
four new edge lengths added as configuration variables. This sector of the 4D dynamics behaves 
very similarly to 3D gravity and was also discussed in [38]. Furthermore, a derivation of the 
symplectic structure for this sector, starting from a first order formulation of canonical discrete 
gravity, may be found in [38]. The results obtained there are in agreement with (|5,6p . The reason 
why this sector is totally constrained is that the 1-4 Pachner moves neither lead to inner triangles 
nor inner edges. Hence, the hypersurface triangulations produced by means of this move allow 
for bulk triangulations without any inner triangles which could carry inner curvature. That is, 
we are considering again just flat triangulations. 

This will, in general, change if in addition we consider 2-3, 3-2, and 4-1 Pachner moves. 
As the 1-4 Pachner move will not lead to any additional inner triangles or edges but adds four 
boundary edges, we can consider it as just adding lapse and shift degrees of freedom. Indeed, 
the four new edge lengths will remain free parameters, at least before doing any other moves. 
The 2-3 Pachner move, on the other hand, generates an inner triangle but no inner edge. We 
shall see that this move can be interpreted as adding a curvature degree of freedom (or physical 
degree of freedom as opposed to the lapse and shift gauge degrees of freedom). The 3-2 move, 
furthermore, produces an inner edge for which the Regge equation of motion — canonically in the 
form of a pre-constraint — must be satisfied. Consequently, this move can be considered as the 
'true evolution step' in the sense that it requires the solution of an equation of motion, involving 
curvature degrees of freedom. Finally, the 4-1 Pachner move generates six inner triangles and 
four inner edges and will feature four pre-constraints. For the same reasons as in the 3D case, 
the pre-constraints will be automatically satisfied if we consider initial data leading to a flat 
solution. For curved solutions, however, these pre-constraints will fix some free parameters 
introduced in the previous steps. 

5.1 The 1—4 Pachner move 

This move is the 4D analogue of the 1-3 Pachner move discussed in section 14.11 Glue a 4- 
simplex a onto a given 3D hypersurface in such a way that its 'bottom tetrahedron' r is 
identified with a tetrahedron in Similarly, consider the situation in which one simplex a 
which shares one tetrahedron r with Ej. is removed from a given triangulation (equivalent to 
adding a tetrahedron with opposite orientation). 

This move acts as a 1-4 Pachner move on E^. It introduces one new vertex and four new 
edges into the new hypersurface, but it does not render any edges internal. The 1-4 Pachner 

9 Once more, we could equally well start by producing an initial hypersurface of more complicated topology. 
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moves replaces the tetrahedron r with the subdivided tetrahedron (consisting of four smaller 
tetrahedra); the boundary of the tetrahedron r does not change and we have the same four 
triangles t e and six edges l e before and after the move. Through the subdivision there will 
appear four new edges e n and six new triangles t n adjacent to these new edges, see figure [TTT 

1-4 





4 1 



Figure 11: The 1-4 Pachner move and its inverse, the 4-1 Pachner move. 

Since there are no edges which become internal, there are no equations of motion to be 
satisfied. As a consequence, the lengths of the four new edges, labeled by n, can be freely chosen. 
That is, one has a fourfold freedom in choosing the 'tip' of a, which can be parametrized by 
lapse N and shift iV a variables. As opposed to the 3D case (and the continuum) this freedom 
will, in general, be restricted by the appearance of pre-constraints in later moves (see also the 
discussion in section [531 about tent moves). 

We proceed in the same way as for the 1-3 move and use the generating function 

Gl-4(4+l>P6 i ^k+l^Pe'i Ik'lk+l) = ^k+lPb + ^ Ik+lPe + S<r{l%+li 'fc+l) > (5-9) 

b e 

where we use the one-simplex action 

5, = £^(^,^(-0^^ .(5.10) 

t e Cv t™C<7 

The equations of motion are then 

% = lk+i , Pt +1 = Pb , (5-11) 

k ~ L k+l i Pe ~ Pe+ o/e ~ Pe ~ 2^ ^ e + 2^ gfe ^ ~ d t n > ^ bAl ' 

fc+1 t e Ccr k+1 t n Ccr k+1 

P k n = V , Pn +1 = ^ = E|^(--^)- 

Ul k + 1 t n Co . Ul k + 1 

As for the 1-3 move the Vj; and remain arbitrary, nevertheless, we can equire = 
Again, we encounter post-constraints for every new edge, namely the second equation in (|5.13p . 
These post-constraints fix the momenta as a specific combination of the exterior angles. At a 
four-valent vertex, which can be seen as the tip of a four-simplex, this combination of exterior 
angles can be expressed as a function of the adjacent length variables (namely the lengths of 
the simplex a). Such a post-constraint will always appear if we produce a boundary edge in the 
course of evolution at which no bulk triangles (as potential carriers of curvature) are hinging. 

In the later discussion we shall see that these constraints are, in general, not preserved by the 
other moves. Additionally, the free data, that is the lenghts of the four new edges, may, in gen- 
eral, become fixed by pre-constraints of later moves. This is related to the fact, that the Regge 
discretization does not preserve the diffeomorphism symmetry of the continuum |1CH [TT] . The 
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constraints and the free data discussed here, therefore, correspond to the Hamiltonian and dif- 
feomorphism constraints and the gauge freedom (of lapse and shift) of the continuum. Through 
the breaking of diffeomorphism symmetry by discretization, the constraints encountered here 
will generally not be preserved by the other Pachner moves but, instead, turned into so-called 
pseudo-constraints \10\ [TTT [T6| I24j . These are equations of motion for canonical data of two 
different times, in which these data are only weakly coupled to each other. (More precisely, the 
eigenvalues of the appropriate Hessian matrix associated to these equations of motion will turn 
out to be small compared to those of Hessians associated to proper evolution equations.) 

5.2 The 2 3 Pachner move 

Next, let us discuss the 2-3 move which turns out to generate a dynamical (or physical) degree 
of freedom. 

Consider the situation in which a 4-simplex a is glued to T,^ in such a way that two of its 
(adjacent) tetrahedra are identified with two (adjacent) tetrahedra in This gluing process 
is only possible when these two tetrahedra in Sj, are not part of the same 4-simplex in the 
underlying triangulation (for the five vertices of these two tetrahedra are already the five vertices 
of a). In contrast to the 2-2 move in 3D, the extrinsic curvature angle around the triangle t° 
along which the two adjacent tetrahedra are identified need a priori (i.e. kinematically) not 
be negative^ On the other hand, we can remove such a simplex from the triangulation (or, 
equivalently, add a simplex with opposite orientation), if the extrinsic angle around t° is positive, 
the two adjacent tetrahedra reside in the same a and there is a piece of triangulation underneath 
Sfc. From the perspective of the hypersurface, this evolution move amounts to a 2-3 Pachner 
move. 




Figure 12: The 2-3 Pachner move and its inverse, the 3-2 Pachner move. 

Two tetrahedra sharing a triangle t° and having a boundary consisting of six triangles t e 
and nine edges e are replaced by three tetrahedra sharing a new edge n and sharing pairwise 
three new triangles t n , see figure [T2J As for all the Pachner moves, the boundary of the three 
new tetrahedra is the same as for the two tetrahedra we started with. Consequently, during 
this move no edge will become internal and, therefore, there will be no additional equations 
of motion which we could solve for the length of the new edge. Furthermore, the triangle t° 
becomes a bulk triangle in the gluing move. This is important, since after the move, the bulk 
triangle can carry a non-vanishing deficit angle, that is curvature. This deficit angle depends 
on the value of the new edge and as this can be freely chosen, we are generating a curvature 

10 Due to the flat embedding in 3D Regge Calculus (without a cosmological constant), the extrinsic angle ip„ 
around the edge e = o in the configuration of the 2-2 Pachner move in section 14.1.11 must be negative for a 2-2 
gluing move to be possible such that eventually e — (for non-degenerate simplices, the dihedral angles are 
always positive). 
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degree of freedom0 This is different from the four free edge lengths that arise in the 1-4 move 
which rather correspond to lapse and shift, and, therefore, gauge variables, in the continuum. 

Because of the new edge at step (k + 1), we extend the phase space at step k by the pair 
(l%,Pn). We use the generating function 

G^ 3 (l b k+1 ,p k b ; lt +1 ,p k e ; = ECiPfc +E l t+iPe + > ( 5 - 14 ) 

b e 

where in this case the one-simplex action is given by 

SaQ>k+iilk+l) = E At<>(lk+iJk+i) + At°(l%+i,lk+i) i 71 ~ 8t°(l%+n ^fc+i)) + 

V A tn (li +1 ,l n k+1 ) (tt - • (5-15) 



t e C(T 
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The equations of motion read 

$ = 4+1 , P 6 fc+1 = , (5-16) 



't°j 



+ E|^(--^), (5-17) 

fc+l t n C(7 + l 

The second equation in (|5.17p determines the momenta p k+l at time k+ 1 as a function of the 
ones at time A;, as well as the new edge lengths l k+1 and = l%- The new lengths l k , l k+1 remain 
undetermined, but we may again choose 1% = Accordingly, we find a post-constraint at 

time (fc + l), namely the second equation in (|5. 18|) . It determines the new momentum variable 
as a function of the edge lengths at step (k + 1). Geometrically this constraint arises as by 
construction there are no bulk, but only boundary triangles hinging at the new edge. Hence, we 
can compute the exterior angle from the boundary geometry alone, i.e. the length variables of 
the hypersurface Sfc+i. 



5.3 The 3—2 Pachner move 

The 3-2 Pachner move is the inverse of the 2-3 move: a 4-simplex a is glued onto T, k in such a 
way that three of its (adjacent) tetrahedra, sharing an edge o, are identified with three adjacent 
tetrahedra in E^. This gluing is possible only when these three tetrahedra reside in three distinct 
4-simplices of the triangulation underlying As in the case of the 2-3 move and, by virtue of 
the possible absence of a flat embedding, the extrinsic angles at the three triangles t° which are 
each shared by two of the three tetrahedra in need a priori not be negative; the values of the 
resulting deficit angles around these triangles will eventually be determined by the dynamics. 
Likewise, one can remove a simplex a from the triangulation underlying (or, equivalently, 

n For instance, it is possible to generate the complete 4D star of a triangle, i.e. rendering it internal, by starting 
out with one 4-simplex and performing a sequence of 1-4 Pachner moves on it and a final 2-3 move. During this 
sequence, the triangle shared by all simplices has become internal, but there are only boundary edges. This star 
of the triangle, therefore, carries curvature without any internal edges. Hence, there are no Regge equations to 
be satisfied. 
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glue a simplex with opposite orientation to E^.), if the three extrinsic angles are positive and the 
three tetrahedra are part of the same 4-simplex. These gluings/removals appear as 3-2 Pachner 
moves, that is three tetrahedra sharing one edge o and pairwise three triangles t° are replaced 
with two tetrahedra sharing one triangle t n , such that the boundary of the two complexes, 
consisting of six triangles t e and nine edges e remains unchanged (see figure fl~3"j) . The edge o will 




Figure 13: The 2-3 Pachner move and its inverse, the 3-2 Pachner move. 

become internal, in consequence we will have to satisfy the equation of motion corresponding 
to this edge. This equation of motion will be implemented via a pre-constraint and momentum 
matching. It may either not be possible at all to satisfy this pre-constraint, or only for specific 
choices of parameters, that is lengths variables, which appeared in previous moves. 

Since we have a variable l° k that does not appear at step (k + 1), we extend the phase space 
at time (k + 1) by the pair (l k+1 , p k+1 ) ■ The generating function for the 3-2 move reads 

G 3 . 2 (l b k ,p k b +l ;li,p k e +1 ;l k ,l° k+1 ) = -J2M +l ~J2 l kPe +1 + S ^ l °k) > (5-19) 

b e 

where the action S a is given by 

J2A t o(l e k ,l° k ) (ir -9^(11,11)) . (5.20) 

t°Ca 

The equations of motion amount to 

l b k+1 = l\ , p\ = p k b +1 , (5.21) 
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(tt-^c) , (5-22) 



P k +1 = , p k Q 
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This time the second equation in (|5.22[) determines the momenta p k+l at step (k + 1) as a 
function of the ones at time k and the old edge lengths I? and It. 

The second equation in (I5.23|) . on the other hand, is a pre-constraint. It requires that the 
momentum p k be given as a function of the edge lengths at time k. This equation is automatically 
satisfied if the edge o has just been created by a 2-3 move, but in general this condition may not 
be satisfied. In this case, one either has to tune the new length variables that arose in previous 
2-3 or 1-4 moves, or perform a different move. 
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5.4 The 4 1 Pachner move 



Finally, let us consider the 4-1 move, which is the inverse of the 1-4 move (see figure [TT]) . A 
complex of four tetrahedra sharing a vertex, four edges o and six triangles t° is replaced by one 
tetrahedron, such that the boundary of four triangles t e and six edges e remains unchanged. The 
4-1 gluing move is only possible if these four tetrahedra are part of four different 4-simplices 
in the underlying triangulation, while the six extrinsic angles at the six t° are kinematically 
unrestricted. From the 4D perspective the four 'old' edges become bulk edges, consequently, 
we will have four equations of motion to satisfy. These take the form of pre-constraints and 
eventually determine the resulting six deficit angles around the t°. The prerequisite for the 4-1 
removal move, on the other hand, is clearly that the four relevant tetrahedra reside in the same 
4-simplex and the six extrinsic angles at the six t° are positive. 

We extend the phase space at time (k + 1) by the four pairs (l^ +l ,Po +1 ) and define the 
generating function 

G A ^{ii P i + \iip^\iix k+1 ) = -EW^-E^+' + ^M) > ( 5 - 24 ) 

b e 

where we use the one-simplex action 

t e C<J t°Cv 

The equations of motion are then 

A = ll+i , Pb = P k b +1 , (5-26) 

l k - i fc+ i , Pe ~ Pe ~ -Qfe- ~ Pe +2^ ~7^ ie ~ 2^ ~7)f~ ^ ~ 0t °> ' ^^'l 

k t e Ccr k t o Ca k 

OL k t n C(J OL k 

The remain undetermined, as these do not appear at all. Again, we can choose 1% = 
to carry these data along as boundary data. The second equation in (|5.28[) represents four pre- 
constraints, that have to be satisfied by the canonical data at time k, so that the 4-1 move can 
be performed. This will, in general, restrict the free parameters which appeared in the previous 
moves. 



5.5 Example: 4D tent moves 

As a specific example, we can also consider tent moves in 4D which, as already explained in 
section 14.21 can be obtained by a succession of Pachner moves. In particular, an A^-valent tent 
move in 4D can be reproduced through a sequence of one 1-4 move, (N — 3) 2-3 moves, (N — 3) 
3-2 moves and one final 4-1 Pachner mover 2 ! 

Here, we will discuss two tent move configurations: the first one at a four-valent vertex leads 
to flat dynamics, while, in contrast to this, the second example involves curvature and the role 
of the pre- and post-constraints differs considerably from that in the first example. 

12 Note that the ordering of the 2-3 and 3-2 Pachner moves is not entirely fixed: one always has to start with 
a 1-4 move to introduce the new vertex v n +i and a subsequent 2-3 move and finish off with a 3-2 move prior to 
the final 4-1 move (which removes v n from E„+i), but the ordering of the 2-3 and 3-2 moves in between can 
be chosen freely (according to the given configuration). 
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To begin with, consider the simplest tent move configuration in 4D, namely the four-valent 
tent move. Starting from a four-valent vertex vq in the boundary surface So of a single 4- 
simplex, perform first one 1-4 and one 2-3 move which introduce five new (freely choosable) 
edge lengths. Finally, perform one 3-2 move to remove one of the edges e(voi), i € 1, ... ,4, and 
one 4-1 move to remove the remaining three edges eiv^j) and the tent pole eiv^vi) from the 
new Si and solve the corresponding five equations of motion or pre-constraints. 

After the 3-2 move the boundary configuration, in fact, corresponds to the configuration 
of a stacked sphere, namely to the configuration of one simplex on which one 1-4 move has 
been performed with the tip of the second simplex pushed "inwards". A stacked sphere is a 
triangulation of the 3-sphere which can be obtained by performing a sequence of 1-4 Pachner 
moves on the 3D boundary surface of a single 4-simplex and therefore necessarily possesses a 4D 
flat interior as there are no internal edges. Any triangulation whose boundary corresponds to a 
stacked sphere configuration must, consequently, be flat as well independent of the existence of 
internal triangles (up to possible discretization artifacts [TU]). 

That is, the 3-2 Pachner move plays a key role here in that its pre-constraint (or equation 
of motion) imposes flatness of the deficit angle around the internal triangle generated in the 
course of the previous 2-3 move and thus establishes one non-trivial condition among the five 
free parameters of the 1-4 and 2-3 moves. Hence, there exists a four parameter family of 
solutions to the equation of motion of the 3-2 Pachner move, all of which correspond to flat 
geometries. Therefore, these solutions also automatically solve the equations of motion of the 
final 4-1 Pachner move (the boundary configuration after the 4-1 move also corresponds to a 
stacked sphere). Thus, the four pre-constraints of the final 4-1 move are automatically satisfied 
and no further non-trivial conditions on the remaining four parameters arise. Accordingly, after 
the 4-1 Pachner move four parameters coordinatizing the position of the vertex remain free and 
one obtains a four-fold gauge symmetry associated to lapse and shift degrees of freedom. 

This situation changes for 7V-valent tent moves with N > 5, since the boundary configura- 
tion at none of the individual Pachner move steps corresponds to a stacked sphere and therefore 
does not necessarily imply flatness. But non-vanishing curvature in Regge Calculus generically 
breaks the vertex displacement gauge symmetry |1CH [XT] and all free parameters become fixed. In 
particular, both the 3-2 and 4-1 move attach internal triangles to the edges in the hypersurface, 
which generally introduces a dependence on data from other evolution steps, e.g. see equation 
(15. 6h . The (post-) constraints after the 1-4 move will, in general, not be preserved by the follow- 
ing Pachner moves. Instead, these transform into pseudo-constraints [UJ. Pseudo-constraints 
can be understood as equations between the phase space data at time k which further depend 
(weakly) on phase space data from previous time steps. 

As the second example consider, for simplicity, a 'symmetry-reduced' tent move at a five- 
valent vertex vq, also used in jlOt 111] , where only two dynamical length variables, a n and b n , 
arise at each tent move step n, apart from the tent pole length. There are six tetrahedra in T,q 
with vertices 

w 124, t> 134, i;o234, v \25, v 135, v 235 . (5.29) 

Accordingly, we will have nine triangles of the form t(vQij) with i,j = 1, . . . 5 in this trian- 
gulation, five edges of the form e(voi) and nine edges of the form e(ij) (all possible ordered 
combinations of i,j € {1, ... 5} with the exception 45). This situation is depicted in figure HH 

'Symmetry-reduction' here means that all the lengths of the boundary edges e(ij) are set to 
1 El and imposing l e ( Vn ^ = a n , i = 1, 2, 3 and Z^™ 4 ) = / e K 5 ) = b n at each tent move step n. 

13 This is possible since the vacuum Regge equations are invariant under global rescalings. 
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Figure 14: Illustration of the symmetry-reduced 3d star of a five-valent vertex v , consisting of the six 
tetrahedra 

The 4-simplices associated to the Pachner moves reproducing this tent move are then all 
of the same type a(voViijn), where i,j take values in 1,2,3 and k in 4,5. Hence, the first 
1-4 move already introduces all three new (freely choosable) parameters a\ , b\ and the length 
of the tent pole between vo and v\ which we will call t\. Specifically, let us consider the 
case where we glue the six 4-simplices onto the tetrahedra (|5.29j) in the following order: 1. 
a{vQV\12A), 2. <r(tioWil34), 3. cr(t>ot>il25), 4. <7(t>o^il35), 5. <r(vofi234), 6. O"(uofi235), which 
corresponds to one 1-4 move, two subsequent 2-3 moves followed by two 3-2 moves and a 
final 4-1 move. These Pachner moves are then to be implemented canonically by equations 
(|5. 1115.1315.1615. 1715.1815.2115. 2215.2315. 26l5.28j) . The necessary action contribution of each of 
these six simplices is of the general form 

S a {v v lijK ) = 2A?(c l 7r-e?) + A b t (c 2 7r-e b t )+A a (c 3 7r-e° a ) 

+A 1 a (c 4 iT-0 1 a ) + 2A o b (c 5 iT-6 o b ) + 2A 1 b (c 6 ir-8 1 b )-A6 , (5.30) 

where 

6®, A® a are the dihedral angle and the area of the triangle A(voij), 
6 b , A® are the dihedral angle and the area of the triangle A(vozk), 
Of, A® are the dihedral angle and the area of the triangle A(voV\i), 
9 b , A\ are the dihedral angle and the area of the triangle A(vqV\k), 
9\, A\ are the dihedral angle and the area of the triangle A(viij), 
6 b , A\ are the dihedral angle and the area of the triangle A{v\in) 
6, A are the dihedral angle and the area of the triangle A.(ijn), 

respectively, and c m G {0, ^, 1}, m = 1, . . . , 6, depending on whether the corresponding triangle 
has been present prior to the move under consideration, is newly introduced or becomes internal 
(see the discussion in the paragraph following (|3.2p in section 13. ip . 

In the sequel we confine our focus to pre- rather than post-constraints since the latter are 
automatically satisfied after the moves. The pre-constraint for e(t>ol) of the first 3-2 move 
coincides with the equation of motion for the ao-edges, the pre-constraint for e(t>o4) of the 
second 3-2 move with that of the fro-edges and the single pre-constraint of the 4-1 move (recall 
the 'symmetry-reduction') is equivalent to the equation of motion for the length t\ of the tent 
pole e(voVi). 
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In the general (not symmetry-reduced) situation, we would have four parameters introduced 
by the first 1-4 move and a further 2 parameters introduced by the two following 2-3 moves. The 
latter two parameters determine the curvature, that is the deficit angles on the triangles which 
become bulk triangles during the completion of these moves. In the symmetry-reduced situation, 
on the other hand, the 1-4 move already introduces all three parameters that are allowed by 
our choice of symmetry-reduction. One of these corresponds to a lapse degree of freedom 
(determining the height of the tent pole), the other two can be interpreted as determining the 
value of the curvature. As we shall discuss, in general all parameters introduced by the 1-4 and 
2-3 moves will be fixed by the pre-constraints of the subsequent 3-2 and 1-4 moves. 

At the level of the momenta, the above symmetry-reduction only holds when a tent move 
step is complete, however, not for the individual Pachner moves into which the tent move can 
be decomposed, since the momenta at the different edges (of identical length) get updated in 
different order depending on the order of the Pachner moves. For this reason, we consider the 
momenta of each of the edges individually, which we denote by p%p® K for the initial data and 
p\-„p\ for those associated to the new edges. (Note that the upper index here counts tent rather 
than Pachner moves.) 

The whole point of elaborating on this example is to demonstrate that, in contrast to the 
four-valent tent move, it is possible to solve only a subset of the equations of motion (or equiv- 
alently, pre-constraints) for ao,bo,ti and violate the remaining ones. This implies that the 
pre-constraints of the five-valent tent move are, in general, independent, not automatically 
satisfied and fix parameters which were previously free. 

For instance, consider the following three examples based on the same initial data ao = 
1.1690, = p° 2 = pi = 1.5088, b = 1.1436, p° 4 = pi = 1.2272: a) tableCQshows the example of a 
situation where only the pre-constraint for ao is solved but not the ones for 6o and ti, b) table [2] 
provides an example where the pre-constraints for ao and 6o are solved, but not the one for t±, 
and c) table [3] demonstrates an example where all three pre-constraints are satisfied. Vanishing 
momenta indicate that a constraint (or equation of motion) is satisfied. In these three examples 
we choose to use the pre-constraints for ao,6o to fix the new lengths ai,£>i, respectively, and 
the pre-constraint for t\ to fix ii0 but note that other choices of which length to fix by which 
constraint are possible. 

In conclusion, after the two 3-2 moves two of the three new parameters (i.e. lengths) are, in 
general, fixed and one has a one parameter family of solutions to the two equations of motion 
of these two moves (with table [2] giving one particular example). However, generally these 
solutions are not automatically solutions to the equation of motion of the final 4-1 move. That 
is, here the three pre-constraints are, in general, independent and completely fix the three free 
parameters of the initial 1-4 move in which case solutions to all moves exist, but not a parameter 
family of solutions and therefore no gauge symmetry arises. Nevertheless, the subset of the one 
parameter family of solutions to the two equations of motion of the 3-2 moves corresponding to 
flat configurations does automatically solve the pre-constraint of the 4-1 move and no further 
non-trivial condition arises on the parameters! 15 ! That is, in this special case one parameter 
remains unrestricted after the 4-1 move (the length of the tent pole) and one obtains gauge 
symmetry. 

14 In fact, prior to the 3-2 moves and the implementation of the pre-constraints the new lengths can, in principle, 
be of any value and only get tuned after imposing the constraints. However, here we set the lengths of the new 
edges to the value determined by one or more pre-constraints already from the start. 

15 There are two independent deficit angles in this symmetry-reduced setup, such that flatness imposes two 
independent conditions among the three free parameters Oi,6i,ti which allow to write, e.g. ai and b\ as a 
function of ti. Due to flatness, by (|5.2[> . the three equations of motion (or pre-constraints) for ao,&o and ti are 
automatically satisfied. 
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Table 1: Momentum updating for an example of the five-valent symmetry-reduced tent move decomposed 
into Pachner moves where only the pre-constraint for ao is eventually satisfied. The rows provide the 
updated momentum values after the Pachner move given in the left column (following the sequence of 
the main text). Initial data as given in the text. Two of the three new lengths introduced in the first 1-4 
move have been chosen as ti = 0.2600, b\ = 1.3400. We impose the pre-constraint (|5.23[) for edge e(i/ol) 
(with length ao) of the first 3-2 move (which translates into p\ = after the move) in order to solve for 
ai, given the other lengths and momenta, and use this value from the start. In this case, a\ = 1.4052. 
(Due to the symmetry, the pre-constraints (|5.28[) associated to e(v 2) and e(«o3) of the 4-1 move are 
then also automatically satisfied and translate into "p\ = = p° after the final move.) On the other 
hand, the pre-constraint for the edge e(i>o4) (with length bo) arising in the second 3-2 move and the 
pre-constraints arising in the 4-1 move for edge e(«o5) (with length bo) and the tent pole (with length 
t\) are violated. 



move 




Pi 


Ps 


Pi 




Pti 


Pi 


Pi 


Pk 


Pi 


p4 


1-4 


2.5119 


2.5119 


1.5088 


0.2082 


1.2272 


5.4351 


-1.0316 


-1.0316 





1.1903 





2-3 


-0.1145 


2.5119 


2.5119 


-2.2703 


1.2272 


3.1184 


1.5135 


-1.0316 


-1.0316 


3.5106 





2-3 


-2.7011 


-0.0747 


2.5119 


-2.2703 


0.2082 


-1.0557 


3.9286 


1.3834 


-1.0316 


3.5106 


1.1903 


3-2 





-0.0747 


-0.0747 


-2.2703 


-2.2703 


-3.3723 


1.3463 


1.3834 


1.3834 


3.5106 


3.5106 


3-2 





-2.7011 


-2.7011 


-1.4644 


-2.2703 


-7.5464 


1.3463 


3.9286 


3.9286 


2.7552 


3.5106 


4-1 











-1.4644 


-1.4644 


-2.1113 


1.3463 


1.3463 


1.3463 


2.7552 


2.7552 



Table 2: Momentum updating for an example of the five-valent symmetry-reduced tent move decomposed 
into Pachner moves where the pre-constraints for both ao, bo are eventually satisfied. Further explanation 
and initial data as given in the caption of table Q] and the text, respectively. In this example, the length 
of the tent pole was chosen as t\ = 0.2000. The pre-constraint for edge e(«ol) (with length ao) of the 
first 3-2 move and the pre-constraint for edge e(i>o4) (with length bo) of the second 3-2 move were 
imposed in order to fix the lengths a\ and b±, given the remaining lengths and momenta. In this case, 
ai = 1.3448, &i = 1.2985. The pre-constraint for the tent pole arising in the final 4-1 move is "almost 
solved", however, still violated (p tl = 0.0003 > 0, while the other constraints are satisfied to order 10~ 12 ). 
The latter is a consequence of the fact that the present configuration yields a near-flat geometry in which 
the vertex displacement symmetry of the flat regime is almost preserved |10l 111] . 



move 




Pi 


Ps 


Pi 




Pti 


Pi 


P L 2 


Pk 


Pi 


Pk 


1-4 


2.1317 


2.1317 


1.5088 


0.7195 


1.2272 


5.1219 


-0.6116 


-0.6116 





0.6253 





2-3 


-0.1145 


2.1317 


2.1317 


-1.3172 


1.2272 


3.1121 


1.5659 


-0.6116 


-0.6116 


2.5562 





2-3 


-2.3208 


-0.0747 


2.1317 


-1.3172 


0.7195 


0.0001 


3.6559 


1.4784 


-0.6116 


2.5562 


0.6253 


3-2 





-0.0747 


-0.0747 


-1.3172 


-1.3172 


-2.0097 


1.4397 


1.4784 


1.4784 


2.5562 


2.5562 


3-2 





-2.3208 


-2.3208 





-1.3172 


-5.1217 


1.4397 


3.6559 


3.6559 


1.3084 


2.5562 


4-1 

















0.0003 


1.4397 


1.4397 


1.4397 


1.3084 


1.3084 



This may be further understood from a covariant point of view, namely an analysis of 
the Hessian of the action (second partial derivatives of the action with respect to t±, a±, b±, t^) 
associated to two tent move steps, which was explored in |10j . There it was shown numerically 
that the Hessian possesses a non-vanishing eigenvalue (corresponding to the lapse degree of 
freedom) in the presence of curvature which approaches zero as the configuration approaches 
flatness. Furthermore, the momenta associated to the edges are in general given by 

_ dS n+ i 

Pe~ , (5.31) 

where S n +i is the effective action of tent move step n + 1 with the length variable of the tent 
pole integrated out. In [10\ [TTj it was shown that the matrix of derivatives of these momenta 
with respect to the l^+i is non-invertible too, if the Hessian is degenerate. In this flat case the 
conditions of the implicit function theorem are violated and we cannot solve for the lengths l^+i 
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Table 3: Momentum updating for an example of the five-valent symmetry-reduced tent move decom- 
posed into Pachner moves where all three pre-constraints of ao,bo,ti are eventually satisfied. Further 
explanation and initial data as given in the caption of table [1] and the text, respectively. The three 
pre-constraints are numerically solved by t\ = 0.3039, ai = 1.4387, bi = 1.3832. 



move 






Pi 


Pi 


Pt 


Pti 


P[ 


Pk 


Pk 


Pi 


Pk 


1-4 


2.2970 


2.2970 


1.5088 


0.7969 


1.2272 


5.4635 


-0.7724 


-0.7724 





0.6121 





2-3 


-0.1145 


2.2970 


2.2970 


-1.3947 


1.2272 


3.3397 


1.5390 


-0.7724 


-0.7724 


2.6485 





2-3 


-2.4861 


-0.0747 


2.2970 


-1.3947 


0.7969 





3.7361 


1.4247 


-0.7724 


2.6485 


0.6121 


3-2 





-0.0747 


-0.0747 


-1.3947 


-1.3947 


-2.1238 


1.4027 


1.4247 


1.4247 


2.6485 


2.6485 


3-2 





-2.4861 


-2.4861 





-1.3947 


-5.4635 


1.4027 


3.7361 


3.7361 


1.3528 


2.6485 


4-1 




















1.4027 


1.4027 


1.4027 


1.3528 


1.3528 



as a function of and . That is, the data at step n does not determine the data at step n + 1 
and, in particular, the free parameters of the 1-4 and 2-3 moves are not all fixed which reflects 
the gauge symmetry. Conversely, when the Hessian is invertible in the presence of curvature, it 
is possible to solve for the l^+i as a function of the lengths and momenta at step n. For the 
symmetry-reduced five-valent tent move discussed here the latter means that, given the aim to 
construct a tent move from the data present at step n, the three a "priori free parameters of 
the 1-4 move are eventually completely fixed by the three pre-constraints of the 3-2 and 4-1 
Pachner moves. There is therefore no gauge symmetry left in the presence of curvature. 

Consequently, in 4D for non-flat configurations, the pre-constraints provide conditions among 
the canonical data (for instance, on the 'free parameters' introduced in the previous 1-4 and 2-3 
moves) and actually implement the equations of motions. In 3D the post-constraints arising 
after the 1-3 moves are preserved by the other moves and exactly match the pre-constraints 
which come with the 3-1 moves. Hence, in 3D we are in the situation where the pre-constraints 
are automatically satisfied by canonical data generated by previous Pachner moves, that is the 
pre-constraints coincide with the post-constraints. 

The situation is different in 4D. Here the post-constraints from the 1-4 and 2-3 moves 
generally do not match up with the pre-constraints of the 4-1 and 3-2 moves. Thus, the 
parameters introduced by the 1-4 and 2-3 moves will, in general, become fixed by the pre- 
constraints of the 4-1 and 3-2 moves. (Of course, whether all parameters get fixed, depends, 
for example, on the number of different moves one is performing.) For initial data leading 
to flat configurations, however, we have seen that the pre-constraints for the 4-1 moves are 
automatically satisfied by the canonical data generated by the previous Pachner moves. This 
converts the four parameters introduced by the 1-4 move to four gauge degrees of freedom (in 
the sense that the solution for given boundary conditions on an initial and final time slice is not 
uniquely determined). 

6 Symplectic structures for discrete canonical dynamics 

The Pachner move dynamics involves phase spaces of different dimensions. Hence, the question 
arises in which sense these maps can be canonical, i.e. preserve a symplectic structure. We will 
discuss this issue in this section by first considering regular systems, and subsequently general 
singular systems, where we have to deal with constraints. In particular, we will introduce in 
more detail the concepts of pre- and post-constraints. Finally, we will consider the situation 
for Pachner moves, and, more generally, for dynamics (generated by an action) which involves 
a change in phase space variables. 
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6.1 Regular systems 

Here we will review the general Lagrangian and Hamiltonian dynamics of discrete systems; for 
regular systems, which we will define shortly, see the exposition [22]. However, we need to 
specifically discuss singular systems, which appear in the Pachner move dynamics (see also the 
discussions in [16J on constraints in discrete dynamics). 

To start with, we will consider regular systems, where by definition constraints will not occur. 
Consider a system with configuration spaces Qk — Q of equal dimension at every time step k. 
We coordinatize Qk by x\ where i is in some index set determined by the dimension of Q. Time 
evolution is generated by the discrete action contributions Sj. = Sk(xi~-i,Xk). These discrete 
actions may be 'effective ones', i.e. arising from Hamilton's principal function and summarizing 
several basic time steps into one effective time step (as, e.g., in the tent moves). 

Consider three consecutive time steps k = 0, 1,2 and the boundary value problem between 
times k = and k = 2. That is, we have boundary data x l and x l 2 and must extremize 

S := S±(x ,xi) + S 2 (xi,x 2 ) (6.1) 

with respect to x\. This gives the equations of motion 

= dSi + dS2 (62) 

dx\ dxi ' 

which we assume to be solvable uniquely for function of xq,x 2 ■ In case that 

det 7T^- + ( 6 - 3 ) 
ax\ 0x2 

we can invert these solutions for (xi,x 2 ) to obtain the Lagrangian time evolution map 

£1 : (x ,xi) H4 (x 1 ,x 2 ) (6.4) 



from Qo x Si to <2i x Q 2 . Condition (I6.3P has to hold for regular systems. 

From the discrete action we can define the so-called Lagrangian one- and two-forms. For 
regular systems these will turn out to be pull-backs via the Legendre transform of the canonical 
one- and two-forms. The significance of the Lagrangian two form is, that it is preserved under 
the time evolution map f)6.4f) . This is the reason for the Hamiltonian time evolution map (in 
phase space) being symplectic with respect to the canonical two-form. 

To start with, we define two canonical one-forms on Q\ x Q 2 from the variation of the action, 

2 {xi,x 2 ) = -^idx\ 



dx\ 
dS 2 

2 (x 1 ,x 2 ) = -Q~r dx 2 » ( 6 - 5 ) 

where we sum over repeated indices Note that dS 2 = 6 2 — 6 2 and, since d o dS 2 = 0, we 
can define the (single) Lagrange two-form 

d 2 S 

n 2 ( Xl ,x 2 ) = d6+ = d0 2 = -^rdx\ Adx> 2 . (6.6) 

dx\dx 2 

This time evolution map is preserved under the time evolution map (|6.4p . i.e. 

fi! = C{n 2 , (6.7) 
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where C\ denotes the pull-back of C\. To see this, consider S from (|6.1|) as a function on Qq x Q\ 
by using for X2 the solutions Xi{xq, x\) of the time evolution map (|6.4p . Taking the exterior 
derivative of S on Q x Qi, we will obtain only boundary terms because the equations of motion 
with respect to the inner variable x\ hold. That is, 

dS(xnXi) ~ ^dx 1 + ^(^dx i + ^dx*\ 

= -0r + C\6t . (6.8) 

Again, as d o (f = and exterior derivatives commute with pull-backs, we also find 

fil = £[fi a • (6.9) 

This argument can be easily generalized to any time step difference (fci,&2)- 

In order to discuss the dynamics in phase space, we introduce the Legendre transformations 

¥ + S k : (x k -i,x k ) (->• (x k , + p k ) = (^k,*T^ 

( dS k 

F~S k : (x k -i,x k ) (x k -i, ~p k ~ 1 ) = [x k -i,-- — -) (6.10) 



dx k -i 

from Qfc-i x Qk to the phase spaces V k -i, V k (these are cotangent spaces of Q k , Q k +i, respec- 
tively). The regularity condition (|6,3p ensures that the Legendre transformations are (locally) 
invertible. The Hamiltonian time evolution 

Hi : (X!,-Pl) ^ {X2, + P2) (6.11) 

can be defined via the Legendre transform from the Lagrangian one fj6.4j) and coincides with the 
time evolution map generated by S2 via the equations (|2.4p . Furthermore, it is straightforward 
to check that the Lagrangian one- and two-forms, (16.51) and (|6.6p . arise from pulling back the 
canonical one- and two-forms, pidx % and dx l Adpi, respectively, with the Legendre transformation 
(|6.10p . This can be used to show that the Hamiltonian time evolution is symplectic, i.e. preserves 
the canonical two form. 



Diagrammetically, 



Qk x Qk+i 




6.2 Singular systems 

Let us now consider the situation where the regularity condition (|6.3p is violated, that is where 
an equal number of left and right null vectors L\ and R\ occurs, satisfying 

L i,J^2 =Q ; J^*4 =0 (6.12) 
dx\dxr 2 dx\dx 2 

in some open neighborhood in Q\ x Q 2 . (To avoid excessively many indices, we will not introduce 
another index numbering the null vectors.) 
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Firstly, we will discuss the consequences for the Lagrangian time evolution which is defined 
as the space of solutions to (|6.2[) , 



= 



dSt dS 2 



(6.13) 



dx\ dx\ 



namely the submanifold in Qq x Qi x Q2 defined by (|6.13p . Given a particular solution, i.e. a 
configuration (xq, x±, X2) satisfying (|6.13p . note that also an infinitesimally displaced configura- 
tion (xq, xi,X2 + ER2) is a solution. Hence, the solution X2 as a function of Xq, x\ is not uniquely 
determined, we rather face a certain amount of arbitrariness. We may call the directions R2 
'gauge directions'. Note, however, that these directions might get fixed in the later time evolu- 
tion steps, that is 'gauge' refers only to the dynamics of this one time step from k = 1 to k = 2. 
In fact, the null vectors L\ and R 2 of the Lagrangian two-form do not necessarily extend to null 
vectors of the Hessian of the action. 

We can still attempt to define a Lagrangian time evolution map from Qq x Qi to Q\ x Q2. 
However, we just discussed that X2 is not uniquely determined, therefore, we either have to fix 
./V gauge parameters, if there are N independent null vectors R 2 , or map to gauge equivalence 
classes. In any case, instead of having the time evolution map onto a 2Q dimensional space, where 
Q is the dimension of configuration space Q, we map at most onto a (2Q — N) -dimensional one. 
This means that either the time evolution map is only defined on some (constraint) submanifold 
of Qo x Qi or that the map is not injective. Also a combination of the two could occur. 

Assume that constraints do not occur, so that we can define a Lagrangian time evolution 
map L\ on the full configuration space Qo x Si- To this end just fix some 'gauge parameters' in 
Qi x 0,2 to determine X2 uniquely as a function of xo and x\. Notice that at this step S(xq, x\), 
i.e. the action evaluated on the solution as a function of the initial data xq, x\, may in some cases 
even depend on this 'gauge' choice. The arguments (I6.8P to (|6.9p showing that the Lagrangian 
two-form is preserved under time evolution still work in the same way as for regular systems. 

However, the Lagrangian two-form will be degenerate, that is possess 2N null directions L\ 
and as i s obvious from the coordinate expression (I6.6P 



Let us turn to the Hamiltonian dynamics. Due to the iV left and N right null vectors (|6.12p . 
the rank of both Legendre transformations 



is given by 2Q — N. Hence, the Legendre transformations are not onto. Rather, the image will 
be given by (2Q — A)-dimensional submanifolds in the two phase spaces V\ and V2 which we will 
call C± and C 2 , respectively. C± is governed by the pre-constraints which need be satisfied on 
V\ so that a Hamiltonian evolution can be defined and is described by the post-constraints 
that are automatically satisfied by the momenta + p 2 after a time evolution from time k = 1 to 
time k = 2 has taken place. 

The N pre-constraints C-f are defined through the equation 



a 2 



d 2 s 2 



dx\ A dx\ 



(6.14) 



dx\dx 3 2 




(6.15) 



= V) 



(6.16) 
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which has to hold for arbitrary x\ and X2- Differentiating equation (|6.16p with respect to X2 
and with respect to x\, respectively, we find the relations 

= (6.1T) 



d~Pj dx\dx\ 
dC{ dC{ d 2 S 2 



dx\ d pj dx{dx\ 



(6.18) 



Since there are ./V independent left null vectors L\ of ^ f 2 we can conclude from (16.171) that 

1 ax^ax^ 

Pin- 

^T = E^(*i,Pi)£i , (6-19) 

where we sum over all TV left null vectors L\ and 7^ (xi,pi) are appropriately chosen coefficients 
(which also depend on the constraint under consideration). Similar relations can be found for 
the partial derivatives of ■ This specifies the gradients of the constraints; the orthogonal 
subspace is then tangential to the constraint hypersurface C[~, respectively ■ 

The Hamiltonian time evolution map %\ will only be defined on the submanifold C-j~. As 
before it can be generated by the discrete action, i.e. 

P = -75— , ^P = -7^— ■ 6.20 

0x1 0x2 

This makes the appearance of pre- and post-constraints explicit. In particular, we can define 
pre- and post-constraints independently from the Legendre transformation. Given a Hamil- 
tonian time evolution map, we declare the pre-image of this map as the submanifold of pre- 
constraints and the image of this time evolution map determines the post-constraints. In the 
previous sections we have given the Hamiltonian time evolution maps for all the Pachner moves 
(and, hence, for all evolution maps involving a sequence of Pachner moves) based on momentum 
updating between the respective constraint surfaces which are defined in this manner. For these 
cases the pre-image and image of the maps could always be specified (but see also section 16.31 
below). 

Furthermore, if we find a solution to the first of the equations (|6.20p for X2 as a function of 
x\, ~~p l , then the infinitesimally displaced configuration X2 + SR2 is also a solution. That is, as 
in the Lagrangian picture, X2 is not uniquely defined. (Note that + p 2 may also change under 
the transformation generated by eR 2 .) Rather, we have N (preliminary) gauge directions i?2- 
The same can be said for expressing x\ as a function of X2, + p 2 where now the 'gauge' directions 
are given by L\. Concretely, we have the gauge displacements 

6 L x\ = L\, S LP j = -—^jL{ . (6.21) 

Note that these vectors are orthogonal to the gradients of the constraints C± , hence, the 'gauge' 
displacements are tangential to the constraint hypersurface Cj~. 

As regards the preservation of the symplectic structure, consider the canonical two-forms 

uji = dx{ A dp) , oj 2 = dx{ A dp 2 . (6.22) 

We can pull these two-forms back via the embeddings i\ : — > V\ and 1^ '■ — > V2 
to two-forms on the constraint surfaces C-f and C^", respectively. To this end, it is conve- 
nient to introduce coordinates {w^ , yjj? , z k o] with k = 1,2 on and C^", respectively. Here 
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the index a = 1, . . . N indicates gauge directions and f3 = 1, . . . , Q — N physical (gauge in- 
dependent) coordinates. These coordinates can be chosen such that for the embedding map 
z.^f : (wi,yi,z ) ' — y (xi,p ) we have 

dx\ _ j dp] _ d 2 S 2 j 

dwf 1 « a ' dw? dx\dx{ 1,a 

M = A r {/3 t 9 Pi - d2S \ N { 
dyf ' <9yf dx\dx{ 

|| = , , (6.23) 

and for i% : (w 2 ,y 2 ,z 2 ) (->■ (x 2 ,p 2 ) 



I. 

8w% ~ z > a ' <M 3x^4 



dy% 3 dyP d4dx J 2 2,/3 

^ - ^ - M j &■■ (6 24) 

Here we introduced the label a to number the null vectors L\ a and .R 2 a - Furthermore, ^, M 2 ^ 

are vectors satisfying „L\ a = and M 2 ^i? 2 Q = for all a, f3. 
The pull-backs of the two forms can then be written as 

= ( L {a dw \ + N U d Vl) A (-^^(Ll al dw?' +N^ l dy()+Nl^ lj dz^ 

j yax 2 ax 2 j 

(6.25) 

The terms with the second derivative of the action S 2 vanish because these are contracted with 
an antisymmetric form. Moreover, the contractions of, for instance, L\ with N\ vanish so that 
we are left with 

(t+)*w 2 = ^Mlpdy^ A M^Sijdz^' . (6.26) 
3 

From these expressions we see that the 'gauge' vetors &L a in (|6.2ip are degenerate directions 
of the two-form (tj~)*u;i, whereas the corresponding Sji a are degenerate directions of (z.^)*^. 
Finally, expressing the Hamiltonian time evolution (|6.20p directly in the coordinates yjr and z^, 
one can check that the Hamiltonian time evolution preserves the pull-backs of the canonical 
forms 

{l-)*uj 1 = 'H* 1 {l+)*uj 2 . (6.27) 
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In this sense Hamiltonian time evolution is a symplectic map. 
Hence, diagrammetically (Ct = Cj^ by momentum matching) 



> Qk x Qk+i 




> C k+l C V k+l ■ 

6.3 Dynamical systems with changing phase space dimensions 

Such singular systems as discussed above appear in the dynamics generated by the Pachner 
moves. There we had to deal with the problem that the number of variables may change 
from one time step to the next. Let us consider a simple example which will highlight the 
general principle. Consider again three consecutive time steps with variables xq,xi,X2- But 
now assume that among the variables X2 there is a 'new variable' xJJ, in particular, that the 
number of variables at time step k = 2 is Q + 1, whereas it is Q at times k = 0,1. Although this 
problem might be well posed as a boundary value problem keeping the data at k = 0, 2 fixed, 
it will, in general, not be possible to turn it into a well posed initial value problem with initial 
data at times k = 0, 1. 

To describe the situation nevertheless by an initial value problem, we extend the configuration 
space at times k = 0, 1 by the configuration variables Xq and x™, respectively. We can interpret 
x~2 as an initial datum which becomes only relevant at time k = 2. To put it into the usual 
initial value problem form, we extend the configuration spaces at earlier times in order to allow 
to already specify the initial datum x?J = x™ = Xq at time k = 0. 

Now the actions Si(xo,xi) and S^O^i)^) will not depend on the variables Xq and x™ as 
these were just introduced for book keeping purposes. We are, therefore, in the situation of 
irregular dynamics both for the time steps between k = and k = 1 and from k = 1 to k = 2. 
Let us assume for simplicity that the dynamics is otherwise regular. 

For the time step from k = to k = 1 the extension does not interfere at all with the dynamics 
of the other variables. We can thus also extend the phase spaces at times k = 0, k = 1 by the 
pairs (Xfc,Pn), however, the Legendre transformations will map onto the constraint hypersurfaces 
Cq = = and = + p\ = 0, respectively. From the previous discussion we know that the 
canonical two-forms restricted to the constraint hypersurfaces will have degenerate directions, 
corresponding to the variables Xq and x", respectively. Also these are 'gauge directions' in the 
sense of the previous discussion. Hence, we can fix these to some arbitrary value, in particular, 
to the one presribed by x?? ■ Note that the 'reduced phase space', i.e. the constraint hypersurface 
modulo the 'gauge direction' (which is the null direction of the restricted canonical two-form) 
coincides with the unextended phase space we started with. 

Next, let us consider the time step from time k = 1 to time k = 2 described on the extended 
phase space V[ and the phase space Vi (which already includes the pair (x2,Pn))- The matrix 

(6.28) 



dx\dx\ 



has (due to our assumption only) one left null vector L\ = b % n . Consequently, there is also a right 
null vector R\ and the Lagrangian two-form (|6.14p possesses these two degenerate directions. 
The Legendre transformations (|6.15p will again not be surjective; for time k = 1 it will rather 
map onto the constraint hypersurface = ^p 1 = 0. This pre-constraint coincides with the 
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post-constraint + p\ = from the previous time step. Moreover, we will have a post-constraint 
C 2 = on the canonical data at k = 2. 

As we have seen, we will also find 'gauge directions' associated to the constraints. At time 
k = 1 this 'gauge' direction (|6.2ip is just given by the coordinate x 2 , i.e. 5lx\ = 5 l n and 5lp} = 0. 
At time k = 2, on the other hand, we have 

8 R x\ = R 2 , 6 R p 2 i = 1 ^ 7 Ri . (6.29) 

ox 2 ox J 2 

The canonical two-form of the extended phase space V[ restricted to the constraint hy- 
persurface p\ = possesses one degenerate direction corresponding to the 'gauge variable' x\. 
Thus, the reduced phase space together with its canonical form coincides with the unextended 
phase space in conjunction with its canonical form. As previously discussed, the Hamiltonian 
time evolution preserves the canonical two-forms restricted to the constraint hypersurfaces. For 
time k = 1 this is just the canonical two-form of the unextended phase space. For k = 2 this 
is the canonical two-form on the phase space including the pair (x 2 ,pfy but now restricted to 
the constraint hypersurface C 2 = 0. It will possess a degenerate direction described by (|6,29p . 
That is, the time evolution embeds the 'smaller phase space' of k = 1 into the bigger phase 
space at k = 2 such that the (restricted) two-forms are preserved. 

The case of 'old variables' x° that disappear at time k = 2 can be discussed in the same way 
by reversing time direction. 

As an extreme example, consider the cases described in sections U] and [5l namely the time 
evolution step from a zero-dimensional phase space to the phase space corresponding to the 
boundary of a tetrahedron or 4-simplex, respectively. The initial phase space is zero dimensional 
and the phase space at time k = 1 twelve- or twenty-dimensional, corresponding to the boundary 
data of a tetrahedron, respectively 4-simplex. In this case the canonical two-forms at k = are 
just the zero forms as the initial phase space is zero-dimensional. Indeed, the time evolution 
preserves this property, as in both cases we had six, respectively ten constraints, fixing all the 
momenta as functions of the lengths variables. Moreover, one can check that the canonical 
two-forms restricted to the constraint hypersurfaces are totally degenerate, that is identically 
zero. 

Finally, let us consider the Pachner moves explicitly. These are rather directly defined in 
the Hamiltonian picture. The reason is that the Lagrangian picture requires the knowledge of 
configurations and velocities at a given time which in the discrete case translates to configuration 
data of two consecutive time steps. For a Pachner move (or any other evolution move, like the 
tent move, which does not cover the entire hypersurface), however, the two consecutive time 
steps overlap, so that we do not have a full set of canonical data. In the Hamiltonian picture 
the velocities, or in the discrete the configuration data of the second time step, are replaced by 
the momenta, which are defined at the same time step as the configuration data. That is, we 
just need one time slice to encode the canonical data. 

For the Pachner moves of type M—N with N > M in 3D and 4D we obtain K new edge 
lengths (K = 1 for the 2-3 move and K = N for the 1-N move). We, therefore, extend the 
phase space at time k by K pairs {l^,Pn) which are matched by the K new edge lengths and 
conjugate momenta {^k+i^Pn +l ) a * time k + 1. Additionally, we have pairs (l b ,p b ) which do not 
change at all during this evolution move and variables (l e ,p e ) for which only the momenta are 
updated. The Hamiltonian evolution map for this kind of Pachner moves (see, for instance, 
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(|5,16rf5TT£j) ) is given by the following momentum updating 

A = t + i , P k b +1 = Vl , (6.30) 

aq lie m \ 

L k — L k+1 ' Pe — Pe H oTi , (O.dl) 

UL k+l 

P k n = , Pt +1 = dSa[l l^ +l) ■ (6-32) 

Here S a stands for the action of the glued tetrahedron or 4-simplex and we choose this action 
to be a function of the length variables at time k + 1. 

Equations (|6.32p contain K pre-constraints, namely C k = p!^, and K post-constraints 

_ „fc+l UJ <r\ l k+V l k+l) fa oq\ 

°n — Pn ) [0.66) 

since this equation only involves canonical data from time step k + 1. Note also that the shape 
(I6.33P of the constraints ensures that these have vanishing Poisson brackets among each other, i.e. 
they form and Abelian (sub-) algebra^ (This is consistent with the previous discussion showing 
that to every pre- or post-constraint we can associate a degenerate direction of the symplectic 
forms restricted to the constraint hyper surf aces. This actually means that the constraints are 
first class.) 

Apart from this, as we discussed for each Pachner move, the Hamiltonian evolution map 
(16.301 [6732]) does not specify all variables uniquely and we rather obtain 'gauge directions' both 
at time k and k + 1. Note that these are 'gauge' only with respect to this evolution move 
from step k to k + 1 and, more precisely, with respect to the symplectic forms restricted to the 
constraint hypersurfaces arising in these moves. Subsequent moves might result in a different 
set of (pre-) constraints with respect to which the original 'gauge directions' do not constitute 
degenerate directions of the restricted symplectic form anymore. 

At time k the length variables Z£ remain undetermined, so that we have Slx\ = 5 l n L^ and 
SiPi = (where i stands for any of the indices b, e, n and are K linearly independent vector 
fields). Also at time k + 1 we have K undetermined variables, namely the edge lengths 
which we referred to as free parameters. These correspond to 'gauge directions' 

a2 a 

S R ll +1 = 6iRl +1 , S RP 1+ l = -. g ° nl R n k+l , (6.34) 

where 1 are K linearly independent vector fields. We have seen that the new edge lengths 
^+1 can be interpreted as additional initial data. These can turn either into physical degrees 
of freedom (as for the 2-3 moves where the new edge length determines the deficit angle at the 
bulk triangle generated by this move) or into proper gauge degrees of freedom (as for the 1-4 
move for initial conditions leading to flat solutions). Additionally, the remain undetermined, 
however, in the previous discussions we pointed out that these can be 'gauge fixed' to be equal 

to q +v 

As for the case discussed in the previous section I6T21 the 'gauge directions' are tangential to 
the constraint hypersurfaces defined by the constraints and C^ +1 , respectively. 
We can consider for both time steps the canonical symplectic forms 

u k = dl\ A dp k b + dl% A dp k e + dll A dp k n , u k+1 = dl b k+l A dp k b +1 + dl e k+l A dp k e +1 + dll +1 A dp k n +1 

(6.35) 



6 Dropping the evolution index k, one directly computes {C„, C7„/} = g f n g" n i — a fn'g"n ~ 0- 
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and restrict these to the constraint hypersurfaces in the respective phase spaces. For the (ex- 
tended) phase space V k at time k we can use the canonical embedding i k : {l k ,p k , l k ,p k , ^) ^ 
(^jfc'Pb' ^kiPn = 0) an d the pull-back of the symplectic form gives just the symplectic form 
of the unextended phase space 

(i k )*co k = dl b k Adpl + dl%Adp h e . (6.36) 

For time step k + 1 consider the embedding 



fc fc+l • I'fc+ljPfc I'fe+ljPe ) 'fc+lJ ^ V'fc+ljPfi i L k+l-,Pe > 'fc+ljPfe+l — — J ' 



In this case, the pull-back of the symplectic form is given by 

(L k+1 )*u k+1 = dll +1 Adp k b +1 + dlt +1 Adp k e +1 + dl% +1 A — ° Q{e dl e k+1 . (6.38) 

Now the time evolution equations (|6.30l - l6T3"2~j) express the canonical coordinates at time (k + 1) 
as function of the canonical coordinates at time k and as function of the length variables l k+1 
(as these are not determined by the time evolution map). Using this for the symplectic form 
(|6.38p at time k + 1, we can conclude that the time evolution map T~i k defined by (|6.30f[6732]) 
preserves the restricted symplectic forms: 

nt^k+iT^k+i = di b k Adp k b +di e k Adp k ; + di e k A ° dit +1 + dii +l h — ° di e k+1 

at fc+i OT fc+i at fc+i at fc+i 

= dl b k Adpl + dl%Adp k e . (6.39) 

The Pachner moves of type N-M with N > M can be discussed analogously, one just has 
to reverse the time direction. 

The 2-2 move in 3D is the only Pachner move for which we extended both the initial and 
final phase space because we had both an 'old' and 'new' variable pair, see section 14.1.11 The 
equations of motion were given by 



ll = Ci > P k b +1 = V\ , (6.40) 

e 

fc+1 

an 



f)C lie 10 in n 

ie 7 e fc+1 k , UlDr \ l k+V 'fc' l k+U fa a->\ 

h = <fc+l > Pe = Pe + HTe > ( 6 - 41 ) 



P k n = , pf = g%^tj , ( 6 .42) 

OL k+i 

P k +1 = , p k = _ gMklM±i) . (6 .43) 

OT fc+i 

Here the only pre-constraint is = and the only post-constraint = 0. Note that 

the second equations in both (16.420 and (|6.43p do not constitute constraints as these depend 
on canonical data of both times k and k + 1. Correspondingly, the only gauge directions one 
encounters here arise from the fact that l k and Z£+i remain undetermined by the equations 
of motion. Therefore, these can be gauge fixed to arbitrary values, for instance l k = l k+1 
(although the latter length needs to be determined from the equations of motion) and l k+ i = l k - 
The canonical symplectic form on the extended phase spaces pulled back to the constraint 
hypersurfaces just gives the canonical forms of the unextended phase spaces. Consequently, we 
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see that in this case the dynamics can be reduced to the unextended phase spaces. The only 
configuration variables that change are l k \- > l k+1 - 

In summary, we have seen that we can define suitably preserved symplectic structures for all 
Pachner moves. To this end we introduced extended phase spaces to circumvent the problem of 
changing phase space dimensions. However, we considered the canonical two-forms restricted to 
the constraint hypersurfaces, and these constraint hyper surf aces, upon factoring out the 'gauge 
directions', coincided with, or were submanifolds in the original, unextended phase spaces. 

An alternative construction is to extend the Hamiltonian time evolution maps, so that these 
are defined on the full (extended) phase spaces. For instance, considering the M-N moves with 
N > M, the relevant changes from (l6.30lJ6T32~j) would be 



*<k+l 

That is, we obtain for all variable pairs % = b,e,n 



P k n +1 = d+Qjf- ■ ( 6 - 44 ) 



The time evolution map is defined on the full extended phase space and its image is given by 
the full phase space at time k + 1. We can then impose the pre-constraints p\ = 0. The post- 
constraints are given as the image of the (pre-) constraint hypersurface under the time evolution 
map. Maps of the form (|6.45j) . which only involve momentum updating via a generating function 
preserve the canonical two-form (of the extended phase space) 1£0 The previous discussion showed 
that also the canonical two-form restricted to the constraint hypersurfaces is preserved. 

Such a construction is also possible for the 2-2 move, i.e. the extended dynamics would also 
be of the form (|6.45p . The pre-constraint p k = and the post-cosntraint p k+1 = 0, however, 
have to be imposed independently. The second equation in (|6.44|) then changes into another 
post-constraint prescribing p^ +1 as a function of the six edge lengths l l k+l of the tetrahedron, 
where i = e, re, o. Alternatively, we can solve for the length l k+1 = l k as a function of the 
five edge lengths Vj^ +1 ,l k+1 and the momentum p^ +1 - Also the second equation of (|6.45p turns 
into a (pre-) constraint determining p k Q as a function of the six edge lengths = e,n,o, or, 
alternatively, l k = l k+l as function of the five edge lengths l° k ,l e k and p k Q . That is, we obtain 
two constraints on each of the extended phase spaces Vk and Vu+i- These are second class 
constraints. Compared to the previous treatment, where we had for, say, V\- only p\ = as a 
constraint and as a gauge degree of freedom, we can consider the second constraint between 
12, 1?, l k and Pg as a gauge fixing condition. This shows that also the induced symplectic forms 
on the (gauge fixed and) constraint hypersurfaces are preserved under time evolution. 



7 Interpretation and role of the constraints 

In spite of the fact that those pre- and post-constraints which arise from momentum updat- 
ing form an Abelian (i.e. first class) Poisson-algebra, they do not, in general, generate gauge 
symmetries of the discrete action. Their phase space flows are associated to null vectors of the 
Lagrangian two-form, but these do not, in general, extend to null vectors of the Hessian of the 
action. However, a degenerate Hessian on solutions implies the existence of flat directions at the 

17 Using the extended map Hk defined by (|6.45[) and uik+i = dll :+1 Adp i +1 , one immediately verifies by symmetry 

and antisymmetry (Hk)*Uk+i = dl\ A d [Pi + ^f) = dl\ A dp*i + a f s °, dl\ A dl\ = dl\ A dp'l = u) h . 

\ ° l k I 9l k dl k 



44 



extrema of the action and thereby the presence of gauge symmetries. Rather, these constraints 
manifest non-uniqueness of solutions and reflect the lack of knowledge at a given time step about 
the full solution (e.g. as a consequence of varying phase space dimensions); in Regge calculus 
at some step k there may not exist sufficient information about the triangulation at step k + 1 
(or vice versa) and a hypersurface may be 'forgetful' about the 'past' or 'future'. For example, 
the post-constraints associated to the new edges simply implement the fact that the lengths of 
these edges can be freely varied in accordance with the canonical data at step k, but not that 
this will eventually be a symmetry of the action. 

On the other hand, the gauge symmetries in Regge calculus are given by those variations 
of the edge lengths which leave the action invariant and correspond to the displacement in flat 
directions of vertices in the bulk of the triangulation. These symmetries are generically broken 
in the presence of curvature [10l [UJ H3J Ej . Such a vertex displacement gauge transformation 
is specifed by the (four linearly independent) null-eigenvectors of the Hessian of the piece of 
action corresponding to the star of the vertex in question [lOl [Til EH [39] . The null-eigenvectors 
- if existent — are, thus, associated to a given vertex and define the flat (gauge) directions in 
which the vertex may be displaced without changing the action. In contrast to this, a priori 
each pre- or post-constraint of the individual Pachner moves is associated to a given edge in the 
hypersurface under consideration, rather than any particular vertex. It describes the variation 
of the corresponding single edge length, however, there is no invariant way of expressing this as 
the displacement of the one or the other vertex at the ends of the edge in some embedding space. 
Nonetheless, in the canonical formalism the variation of the edges can be invariantly associated 
to the displacement of a given vertex v, if the constraints of those edges in the hypersurface 
which are adjacent to v are contracted with the four null-eigenvectors associated to v [UJ@ 




this implements a projection in gauge directions and the correct number of gauge generators. 
At this stage, the contracted constraints generate the gauge transformations defined by the 
null-eigenvectors . 

This was discussed for tent moves in [IT], but will be studied in detail and more generally in 
a forthcoming paper [23] , where we will apply the present formalism to the 4D linearized theory, 
based on an expansion around a flat solution to linear order. This linearized regime, in fact, 
inherits the gauge freedom of the flat background solution and is relevant for the continuum 
limit where the diffeomorphism symmetry of classical general relativity ought to be restored and 
geometries are locally flat. This is where null-eigenvectors of the Hessian necessarily arise which, 
in turn, can be employed to construct proper gauge transformation generating constraints which 
are associated to vertices. 

8 Discussion and outlook 

In the present article we have devised a general canonical formalism for discrete systems involving 
varying phase space dimensions. To this end we had to extend the general theory of discrete 
dynamics from the regular to the irregular case in which constraints necessarily arise. The 
basic ingredient of the new formalism is Hamilton's principle function which may be used as a 
generating function (of the first kind) for a canonical time evolution map. This time evolution 
is generally only defined on a constraint surface and maps to another constraint surface; we 
refer to the pre-image of the time evolution map as pre-constraint surface, while its image is 
called post-constraint surface. The pre-constrainsts, on the one hand, constitute the conditions 

18 Edges internal or not yet occurring at a given step do not matter because they come with vanishing momenta 
and solved equations of motion. 
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that need be satisfied in order for the following time evolution step to be allowed. On the other 
hand, the post-constraints are relations among the phase space data that are automatically 
realized after an evolution step, however, may provide non-trivial conditions to be satisfied by 
subsequent time steps. As shown in general, the symplectic forms restricted to the constraint 
surfaces are preserved under such time evolution maps. 

We have applied the new canonical formalism to simplicial gravity, in particular, and thereby 
constructed a general and fully consistent formulation of (Euclidean) canonical Regge calculus. 
However, the formalism is general and, therefore, adaptable to other discretization schemes (with 
additive action) . The general idea is to implement a canonical forward time evolution by gluing 
single simplices step by step to a bulk triangulation, while backward evolution is simply achieved 
by removing single simplices from the triangulation. The prerequisite for this evolution scheme is 
additivity of the action. The (D — l)-dimensional hypersurfaces then evolve in a discrete 'bubble' 
time through the full L>-dimensional Regge solution akin to the evolution of hypersurfaces in 
canonical general relativity. Both the gluing and removal moves allowed in the present scheme 
can be interpreted entirely within the [D — l)-dimensional triangulated hypersurfaces as Pachner 
moves. Pachner moves are an elementary and ergodic class of piecewise-linear homeomorphisms 
which can map between any two triangulations of a given topology by finite sequences (in 
fact, one could also implement 'spatial' topology changes — aka births of baby universes — by 
generalizing the class of allowed gluing and removal moves, see section [3]). 

The central issue to be addressed when implementing the Pachner moves in a canonical for- 
malism is the varying phase space dimension since Pachner moves introduce and remove edges. 
This obstacle can be handled in a straightforward manner by suitable phase space extensions 
which are controlled by constraints that arise naturally in this framework and constitute the 
canonical incarnation of the equations of motion of edges to the 'future' or 'past' of a given hy- 
persurface. In this fashion, phase spaces before and after any Pachner move can be brought to 
equal dimension, the moves can be implemented via momentum updating as canonical transfor- 
mations and the symplectic forms restricted to the constraint surfaces are preserved under these 
moves. In fact, in an alternative construction without immediate use of Hamilton's principal 
function, the elementary time evolution map can also be extended to the full extended phase 
spaces (i.e. disregarding any constraints) via momentum updating and preserves the full sym- 
plectic structure. This requires, in addition, to impose the pre- and post-constraints separately, 
since momentum updating by itself only specifies the difference of the momenta before and after 
an elementary move, however, not the specific values. 

The implementation of the individual Pachner moves follows one and the same recipe, how- 
ever, each of the Pachner moves assumes a special role in the evolution of the triangulation. 
In 3D the 1-3 move introduces one vertex and equips it with three constraints which generate 
vertex displacements in the flat embedding space; the three new edge lengths thereby assume a 
lapse and shift type gauge parameter interpretation. The pre- and post-constraints of the 2-2 
movj^l simply preserve and implement flatness, while the pre-constraints of the 3-1 move are 
automatically satisfied and this move annihilates one vertex and three constraints associated 
with it. The situation in 3D is special in that all vacuum solutions are necessarily flat and the 
3D Regge action is actually a so-called 'perfect action' [371 HO], i-e. despite being a discrete 
action it, nevertheless, reproduces the continuum dynamics and symmetries and is triangulation 
independent. This has significant repercussions for the initial value problem in the discrete: 
whereas usually the phase space corresponds to the space of solutions this cannot, in general, 
be expected when the phase space dimension varies and here sufficient initial data is rather 

19 Recall that these constraints, in contrast to those from all other moves in 3D and 4D, are only constraints on 
the extended phase space. 
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specified in the course of evolution. Notwithstanding this issue, in 3D the space of initial data 
does correspond to the space of solutions (modulo gauge) because there are only global physical 
degrees of freedom (and the Pachner moves do not change the 'spatial' topology), while the bulk 
is necessarily flat and there exists a three-fold gauge symmetry at each vertex in the triangula- 
tion. The latter symmetry corresponds to lattice diffeomorphisms and is a consequence of the 
action being perfect. 

In 4D, however, the situation changes because of local degrees of freedom and generically 
broken (diffeomorphism) symmetries; the 4D Regge action is not perfect and we face the issue of 
non-uniqueness of solutions in the sense that different solutions arising from the same initial data 
cannot be mapped into each other by gauge transformations because of the broken symmetries 
and are thus inequivalent. As an extreme example of this non-uniqueness for geometries with 
spherical 'spatial' topology take the boundary of a 4-simplex as initial slice which fits into any 
such triangulation and whose data thus by no means specifies the ensuing geometry (or, even 
more extremely, consider the 'no boundary proposal' for triangulations). This is also manifested 
in the particular roles of the individual Pachner moves in 4D: the 1-4 movj^l introduces a vertex 
along with four associated post-constraints which are automatically satisfied after the move and 
reflect the a priori free choice of the four new edge lengths. The 2-3 move generates a bulk 
triangle while only adding a new boundary edge which comes with a new post-constraint. Again, 
the length of this new edge, and thus the value of the curvature angle around the internal triangle, 
is a priori free; the move thereby introduces a (free) curvature degree of freedom. These first two 
moves are the only Pachner moves for 4D simplicial gravity which introduce new edges, however, 
at the same time do not invoke any equations of motion and, consequently, the new edge lengths 
are a priori free parameters. In addition, the 3-2 move is only allowed if the pre-constraint 
corresponding to the equation of motion of the edge that is rendered internal in the course of 
the move is satisfied. Likewise, prior to the 4-1 move which annihilates one vertex, the four 
pre-constraints associated to the four edges adjacent to the vertex which become internal must 
be fulfilled. The pre-constraints of these last two moves impose non-trivial conditions which, 
in general, require a tuning of the free parameters of previous moves (see also the example of 
the symmetry-reduced five-valent tent move in 4D in section 15 . 5|> . In short, in 4D all new 
edge lengths are a priori free data which, however, may become fixed a posteriori by the pre- 
constraints of subsequent 3-2 and 4-1 moves. 

As regards the occurrence of gauge symmetries, we have argued in section \7\ that the pre- 
and post-constraints of the Pachner moves generally do not generate gauge transformations. 
However, the vertex displacement gauge symmetry of Regge calculus exists on flat configurations 
and for this regime, the post-constraints of the 1-4 moves — in analogy to the 1-3 moves in 
3D — do constitute the gauge generators. In general, it depends on the particular configuration 
and, in particular, on choices of data during subsequent moves, whether the 1-4 move introduces 
gauge generators. A specific example of this was treated in section 15.51 

The formalism introduced here is very general and allows to also consider evolution maps 
in which arbitrary pieces of D-dimensional triangulations are glued onto (D — 1) -dimensional 
hypersurfaces (under certain regularity conditions). These can be built up by Pachner moves, 
but may also be considered in an 'effective' picture, where one time step would involve a more 
complicated piece of triangulation than just a simplex. For the example of the so-called tent 
moves we have discussed the 'effective evolution' in 3D in section 14.21 and in 4D in section 15.51 

The new general canonical formalism for discrete systems involving varying phase space di- 
mensions, clearly, warrants further research. In particular, it remains to be investigated under 

20 Here we focus solely on gluing moves. 
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which precise conditions a priori free parameters in 4D canonical Regge calculus remain free 
throughout evolution, i.e. are true gauge parameters, or, equivalently, under which conditions 
post-constraints of previous steps match pre-constraints of later steps (but involving the same 
hypersurface) . This is a prerequisite for providing a constraint classification with clear distinc- 
tion between proper equations of motion and (pseudo-) constraints generating (broken) gauge 
symmetries. The latter question stresses the close connection between constraints and equa- 
tions of motion, which often can be transformed into each other by relabeling time steps (see 
the discussion in section [2]) or by changing the slicing into equal time hypersurfaces. Therefore, 
the question arises whether different choices of 'effective' time steps (or slicings of the bulk 
space-time) may influence the constraints associated to a given hypersurface. This question will 
be important for quantization as one usually attempts to find a quantization of the constraint 
hypersurfaces. Furthermore, the new formalism can also be viewed as a specific algorithm to 
generate solutions and might thereby be interesting for numerical investigations. 

In a forthcoming paper [23], we will address some of these issues explicitly in the linearized 
theory which inherits the gauge symmetry of the flat background. The specific role of each of the 
Pachner moves becomes clear in this regime, namely: the 1-4 move generates four gauge degrees 
of freedom and the same number of gauge generating constraints, the 2-3 move generates one 
graviton, the 3-2 move imposes the only proper, non-trivial equation of motion and annihilates 
one graviton and, finally, the 4-1 move annihilates four gauge degrees of freedom and the same 
number of gauge generating constraints which are automatically satisfied prior to the move. 
This allows for an explicit counting of both gauge and observable degrees of freedom at each 
evolution step in a scenario of an evolving 'spatial' hypersurface. For reasons of continuity, the 
role of the individual Pachner moves must also be the same in a neighbourhood of the linearized 
regime. 

In view of Loop Quantum Gravity (LQG) which involves changing spatial graphs and Spin 
Foam models which employ triangulations, it seems promising to attempt a quantization of the 
Pachner moves and implement them as (possibly unitary) transformations on some Hilbert space 
of 'spatial' geometries. This may help in investigating the connection between LQG and Spin 
Foams (see also [41] for recent advances from the Spin Foam side); due to the fact that the 
classical canonical formalism is based on using the action as a generating function, a canonical 
quantization thereof should coincide with a covariant path-integral quantization of the theory. 
Usually, a canonical program proceeds by quantizing the phase space (or rather the constraint 
surface); here the additional challenge emerges of how to incorporate the fact that initial data 
continues to be specified during evolution and certain a priori free parameters become fixed 
later on. These questions are related to the 'general boundary proposal' [12], which attempts to 
generalize the notion of Cauchy hypersurfaces to arbitrary surfaces for quantum field theories. A 
quantization of the Pachner moves may also shed some light on the 'e ±i5 vs. cos <S" debate in the 
Spin Foam community: because gluing and removal moves add and subtract action contributions 
of simplices, respectively, an amplitude with a single exponential of the action should correspond 
either to pure forward or backward evolution, while an amplitude with the cosine of the action 
incorporates a superposition of both directions of evolution. 

A final interesting question arises: can one construct from the ideas presented here a canon- 
ical formulation of Causal Dynamical Triangulations (CDT) (sec also |43j)? Given that a priori 
all new lengths in 4D canonical Regge calculus can be freely chosen and, in particular, all be 
fixed equal to one while the births of baby universes are disallowed, can one construct 'classical 
CDT solutions' from this, such that the pre- and post-constraints determine the connectivity 
by rejecting or accepting certain Pachner moves? Presumably, this will generically not work and 
one might have to allow, instead, for the lengths to take value in a small interval (1 — e, 1 + e). 
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It would be interesting to investigate whether a quantization with such restrictions can be con- 
nected to CDT. This also requires an extension of the present scheme to Lorentzian signature and 
an implementation of the condition that spacelike hypersurfaces remain spacelike throughout 
evolution. 
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